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Report  Summary 


This  final  report  summarizes  activities  and  work  completed  between  May, 
1990  and  May,  1994. 

During  this  period  we  have  participated  in  eleven  scientific  conferences 
and  advanced  workshops  concerning  various  aspects  of  wavelets  and  signal 
processing.  At  each  of  these  meetings  we  have  presented  the  results  of  our 
investigations  under  this  grant.  As  another  output  of  these  investigations 
we  have  produced  23  original  scientific  articles  on  various  aspects  of  wavelets 
and  their  applications.  For  the  record  we  fist  these  articles. 

1.  Multiresolution  analysis,  Haar  Bases,  and  Self- Similar  Tilings  of  , 
IEEE  Trans.  Info.  Theory,  38  (1992),  556-568. 

Grochenig  and  Madych 

2.  Translation  Invariant  Mtdtiresolution  Analyses, 

Recent  Advances  in  Fourier  Analysis  and  Its  Applications,  Byrnes  and 
Byrnes,  eds.,  NATO  ASI  Series  C,  Vol.  315,  Kluwer,  Dordrecht,  1990, 
455-462. 

Madych 

3.  Wilson  Bases  and  Modulation  Spaces, 

Math.  Nachr.,  155  (1992),  7-17. 

Feichtinger,  Grochenig,  and  Walnut 

4.  Non-Orthogonal  Wavelet  and  Gabor  Expansions  and  Group  Represen¬ 
tations, 

Wavelets  and  Applications,  R.  Coifman  et  al.,  eds.,  Jones  and  Bartlett, 
Boston,  1992,  359-398. 

Feichtinger  and  Grochenig 

5.  Image  Reconstruction  in  Hilbert  Space, 

Mathematical  Methods  in  Tomography,  G.T.  Herman,  A.K.  Louis,  F. 
Natterer,  eds..  Lecture  Notes  in  Mathematics  1497,  Springer  Verlag, 
Berhn,  1991,  15-49. 

Madych 

6.  Wavelets  and  Generalized  Box  Sphnes, 

Appl.  Analysis,  44  (1992),  51-76. 

Lorentz  and  Madych 
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7.  Spline  Wavelets  for  Ordinary  Differential  Equations, 

GMD  technical  report  No.  562. 

Lorentz  and  Madych 

8.  Multiresolution  analyses,  Tiles,  and  Scaling  Functions, 

Probabilistic  and  Stochastic  Methods  in  Analysis  with  Applications,  Byrnes 
et  al,  eds.,  NATO  ASI  Series  C,  Vol.  375,  Kluwer,  Dordrecht,  1992, 
233-243. 

Madych 

9.  A  Riesz  Basis  for  Bargman-Fock  Space  Related  to  Sampling  and  Inter¬ 
polation, 

Arkiv.  f.  Math.,  30  (1992),  283-295. 

Grochenig  and  Walnut 

10.  Reconstruction  Algorithms  in  Irregular  Samphng, 

Math.  Comp.,  59,  (1992),  181-194. 

Grochenig 

11.  Efficient  Algorithms  in  Irregular  Samphng  of  Band-Limited  Functions. 
Proc.  Int.  Phoenix  Conf.  on  Computers  and  Communication  1991, 
IEEE  Comp.  Soc,  1991,  490-495. 

Grochenig 

12.  Some  elementaxy  properties  of  multiresolution  analyses  of 
Wavelets  -  A  tutorial  in  Theory  and  Applications,  C.  K.  Chui,  ed.. 
Academic  Press,  Boston,  1992,  259-294. 

Madych 

13.  Self-similar  Lattice  Tihngs, 

J.  Fourier  Analysis  Appl.  to  appear. 

Haas  and  Grochenig 

14.  Miscellaneous  Error  Bounds  for  Multiquadric  and  Related  Interpola¬ 
tors, 

Computers  Math.  Applic.,  Vol.  24,  No. 12  (1992),  121-138. 

Madych 

15.  Gabor  Wavelets  and  the  Heisenberg  group:  Gabor  expansions  and  short 
time  Fourier  transform  from  the  group  theoretical  point  of  view. 
Wavelets  -  A  tutorial  in  Theory  and  Applications,  C.  K.  Chui,  ed.. 
Academic  Press,  Boston,  1992,  359-298. 

Feichtinger  and  Grochenig 
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16.  Sharp  results  on  random  sampling  of  band-limited  functions, 

Probabilistic  and  Stochastic  Methods  in  Analysis  with  Applications,  Byrnes 
et  al,  eds.,  NATO  ASI  Series  C,  Vol.  375,  Kluwer,  Dordrecht,  1992, 
323-335. 

Grochenig 

17.  A  discrete  theory  of  irregulax  sampling, 

Lin.  Alg.  Appl,  193  (1993),  129-150. 

Grochenig 

18.  Irregulax  sampling  of  wavelet  and  short  time  Fourier  transforms, 

Gonstr.  Approx.,  9  (1993),  283-297. 

Grochenig 

19.  Acceleration  of  the  Frame  Algorithm, 

IEEE  Trans,  Signal  Proc.,  41  (1993),  3331-3340. 

Grochenig 

20.  The  recovery  of  irregularly  sampled  band-limited  functions  via  tem¬ 
pered  splines, 

J.  Funct.  Anal.  125  (1994),  201-222. 

Lyubaxskii  and  Madych 

21.  Orthogonality  criteria  for  compactly  supported  scaling  functions, 

Appl.  Comp.  Harm.  Anal.,  1  (1994),  242-245. 

Grochenig 

22.  Orthogonal  Wavelet  Bases  for  L^{R^), 

Fourier  Analysis:  Analytic  and  Geometric  Aspects,  Bray  et  al.,  eds., 
Marcel  Dekker,  New  York,  1994,  243-302. 

Madych 

23.  Scaling  Functions  and  Sequences  Associated  with  Orthonormal  Wavelets, 
Houston  J.  Math,  to  appear. 

Dlin  and  Madych 

We  bring  attention  to  the  fact  that  twenty  of  these  articles  have  already 
appeared  in  various  scholarly  journals  and  books.  The  first  twenty  one  ar¬ 
ticles  on  this  list  together  with  detailed  summaries  have  been  already  sub¬ 
mitted  with  earlier  technical  reports.  The  final  two  article  on  this  list  are 
appended  at  the  end  of  this  report  and  are  summarized  below. 

•  Orthogonal  wavelet  bases  for  L^(i2”).  An  orthogonal  wavelet  basis  con¬ 
sists  of  dilates  and  translates  of  one  function  or,  more  generally,  a  rela¬ 
tively  small  finite  number  of  functions.  In  this  presentation  we  give  the 
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details  of  the  construction  of  such  bases  where  the  dilation  is  defined 
via  a  fairly  general  Hnear  transformation.  A  Hberal  number  of  examples 
are  given  to  illustrate  the  flexibility  of  the  recipe.  The  significance  of 
this  work  lies  in  the  fact  that,  unlike  the  dyadic  case,  in  the  case  of 
a  more  general  dilation  matrix  there  is  no  prescription  for  construct¬ 
ing  wavelets  from  a  given  scaling  function  or  sequence.  In  this  article 
we  introduce  several  methods  for  constructing  wavelets  with  desired 
properties  in  the  case  of  general  dilation  matrices. 

•  Scaling  Functions  and  Sequences  Associated  with  Ortkonormal  Wavelets. 
It  is  well  known  that  the  integer  translates  of  the  scaling  function  as¬ 
sociated  with  a  given  scaling  sequence  may  fail  to  be  mutually  orthog¬ 
onal.  In  this  article,  we  address  several  technical  questions  related 
to  this  phenomenon.  For  instance,  we  show  that  the  scaling  function 
naturally  associated  with  a  finite  scaling  sequence  always  generates  a 
multiresolution  analysis  and  give  elementary  but  non-trivial  examples 
of  scaling  sequences  which  give  rise  to  pathological  scaling  functions. 
The  significance  of  this  work  lies  in  the  fact  that  it  is  very  important 
in  various  applications  to  know  whether  an  apparent  scaling  sequence 
will  generate  orthogonal  wavelets  via  the  usual  paradigm. 

In  the  previous  technical  report  we  mentioned  that  we  were  preparing  a 
paper  on  results  concerning  the  breakdown  of  the  so-called  scaling  functions 
into  more  elementary  building  blocks.  Since  that  report  we  have  discovered 
that  these  building  blocks  are  closely  related  to  the  class  of  distributions 
which  are  invisible  or  undetectable  by  the  corresponding  family  of  wavelets. 
(Namely,  if  /  is  such  a  distribution  and  {i’j^}j,kez  is  the  corresponding  family 
of  orthogonal  wavelets  then  the  scalar  products  zero  for  all  j  and 

k.)  We  believe  that  these  observations  can  be  developed  into  a  significant 
theory  providing  a  deeper  understanding  of  wavelets  and  the  multiresolu¬ 
tion  analysis  paradigm.  The  principal  investigator  intents  to  continue  his 
investigations  into  this  area. 
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Orthogonal  wavelet  bases  for  LP'^RP') 

W.  R.  Madych* 


Abstract 

An  orthogonal  wavelet  basis  consists  of  dilates  and  translates  of 
one  function  or,  more  generally,  a  relatively  small  finite  number  of 
functions.  In  this  presentation  we  give  the  details  of  the  construction 
of  such  bases  where  the  dilation  is  defined  via  a  fairly  general  hnear 
transformation.  A  liberal  number  of  examples  are  given  to  illustrate 
the  flexibility  of  the  recipe. 


1  Introduction 

1,1  Overview 

In  this  talk  we  present  the  notion  of  orthogonal  wavelet  bases  of 
relative  to  a  general  dilation  matrix  A  and  outline  the  natural  scheme  for 
constructing  such  bases.  This  scheme  was  originally  developed  by  S.  Mallat 
and  Y.  Meyer,  see  [31,  33]. 

An  orthogonal  wavelet  basis  for  L^{1R)  relative  to  dyadic  dilation  is  a 
complete  orthonormal  system  of  the  form 

(1) 

where  ^  is  an  appropriate  function  in  L^{1R)  and  the  indicies  k  and  j  run 
through  the  integer  lattice  2Z  =  {0,  ±1,  ±2, . . .}.  The  elements  of  such  a  basis 
are  often  referred  to  as  wavelets  and  the  function  il;  is  sometimes  referred  to 
as  the  fundamental  wavelet. 

‘Department  of  Mathematics,  University  of  Connecticut,  Storrs,  CT  06269.  Partially 
supported  by  a  grant  from  the  Air  Force  Office  of  Scientific  Research,  AFOSR-90-311. 
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Figure  1:  The  sine  wavelet  and  its  Fourier  transform 


1.2  Examples 

1,2.1 


Consider  the  classical  example 


ifO<a;<  1/2 
if  1/2  <  a:  <  1 
otherwise. 


The  corresponding  system  (1)  is  the  well  known  Haar  basis  for 


1.2.2 

Another  example  is  the  case 

sin  27rx  —  sin  nx 
'ip{x)  - - 

^  TTX 

which  is  related  to  the  theory  of  band  limited  functions  and  cardinal  series. 
Since  the  Fourier  transform  of  ip  is  the  characteristic  function  of  the  set 
:  Tf  <  1^1  <  27r}  it  follows  from  routine  functional  analytic  and  Fourier 
transform  techniques  that  the  corresponding  set  (1)  is  a  complete  orthonor¬ 
mal  system  for  L^(iR).  See  Figure  1  . 

Because  of  the  close  relationship  of  ip  to  the  classical  cardinal  sine,  also 
known  as  the  sine  or  sinus  cardinalis, 

sin  7CX 

sinefa:)  = - , 

TZX 


Figure  2:  The  function  ^  and  the  wavelet  ip  of  Example  1.2.3^ 
we  refer  to  this  basis  as  the  sine  wavelet  basis. 

1.2.3 

Finally  consider  the  function  whose  graph  is  the  second  plot  given  in  Figure 
2.  This  function  is  defined  by 

1p{x)  =  S3(p{x)  —  S2(p{x  —  1)  +  Si(p{x  —  2)  —  SQ(p{x  —  3) 

where  <p  is  the  unique  compactly  supported  solution  of 

(p{x)  =  SQ(p{x)  +  Si(p{x  —  1)  +  S2<p{x  —  2)  4-  S3(p{x  —  3) 

with  L^{1R)  norm  one,  see  the  first  plot  in  Figure  2,  and 

,  ,  fl  +  \/3  S  +  VS  3-V3  l-y/3\ 

{S0,Sl,S2,S3}  -  I  ^  7  ^  )  4  5  4  j 

The  function  ip  is  &  member  of  a  class  of  fundamental  wavelets  introduced 
by  Ingrid  Daubechies  [12]. 


1.3  Dilation  matrices  and  wavelets 

In  the  description  (1)  note  the  role  played  by  the  translations  x  x  +  j, 
j  E  Z  and  the  dilation  x  ^  2x  which  is  expansive  and  maps  the  lattice  Z 
into  itself.  Thus,  more  generally,  if  F  is  a  lattice  in  iR”,  and  A  is  an  expansive 
linear  transformation  on  JR"  which  leaves  F  invariant  then  a  wavelet  basis 
associated  to  (F,  A)  is  a  complete  orthonofmal  system  for  L‘^{Sr)  of  the  form 

(2)  {a'‘^%{A'‘x  - 

where  V'l?  •  •  ■  >^m  an  appropriate  collection  of  orthonormal  functions  in 
L^{1R'^),  a  =  |det>l|,  the  index  k  runs  through  the  integers  Z,  and  the 
parameter  7  runs  through  the  lattice  F.  Wavelets  are  simply  the  members 
of  such  a  basis  and  the  collection  is  the  set  of  fundamental 

wavelets. 


1.4  Examples 
1.4.1 


In  the  case  n  =  1,  T  =  Z,  and  >1  =  3  consider  the  functions  V'/,  I  =  1)2, 
defined  by 


a;)  =  I 


uik  if  k/3  <  X  <  {k  +  1)/Z,  A:  =  0,1,2 
0  otherwise 


where  the  w/it’s  are  scalars  chosen  so  that 


_1_ 


/I  1  1  \ 

«10  Mil  ^12 

^  M20  M2I  M22  ) 


is  an  orthogonal  matrix.  This  definition  implies  that  the  system  of  functions 

{3*'/^^/(3^x-i)},e{i,2}  .fcezjez 

is  an  orthonormal  system  in  F^(JR);  that  it  is  complete  follows  from  an 
argument  similar  to  that  used  to  show  the  completeness  of  the  Haar  system. 
Here  Z  =  1, 2  and  both  k  and  j  run  through  Z. 
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More  generally  if  N  is  an  integer  greater  than  one  and 
/I  1  1  \ 


1 

Vn 


«1,0  ^1,1 


V  Un-1,0  Uff-1,1  •  *  *  UN-1,N-1  J 

is  an  orthogonal  matrix  then  the  collection  of  functions  {tpi,..  .,‘tpN-i}  de¬ 
fined  by 


tpiix) 


"10 


k  iik/N  <x<{k  +  l)/N,  k  =  0,1,. . .  ,N  -  1 
otherwise 


for  1  =  1, . . .  ,iV  -  1  is  a  set  of  fundamental  wavelets.  The  corresponding 
wavelet  basis  of  is  the  collection 

where  /  =  1, . . . ,  iV  -  1  and  both  k  and  j  run  through  ZZ. 


1.4.2 

In  the  case  n  =  2,  F  =  2Z^,  and 


consider  the  function 

•0(a;)  =  2  sinc(Aa;)  —  sinc(a:) 


Where 


and  sinc(x)  =  sinc(a:i)  sinc(x2)  . 


The  function  ^  is  a  bi- variate  analog  of  the  fundamental  wavelet  in  Example 
1.2.2.  The  Fourier  transform  of  is  the  characteristic  function  of  the  region 
:  max{|^il,  |6|}  >  tt,  |6|  +  <  27r},  see  Figure  3.  That 


—  i)}fcezjez2 


is  a  wavelet  basis  of  follows  from  reasoning  identical  to  that  used  in 

Example  1.2.2. 
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Figure  3:  The  wavelet  ip  and  the  support  of  ip  of  Example  1.4.2 

1.5  Contents  and  Notation 

As  mentioned  in  the  overview,  a  natural  framework  for  the  construction  of 
a  wavelet  basis  called  a  multiresolution  analysis  was  developed  by  Y.  Meyer 
and  his  collaborators,  see  [20,  24,  31,  32,  33].  In  view  of  form  of  these  bases, 
(2),  it  should  not  be  surprising  that  the  main  elements  of  such  analyses  are 
the  notions  of  dilation  and  translation  relative  to  a  lattice. 

In  Section  2  we  give  an  outline  of  the  theory  and  include  several  examples 
of  mutiresolution  analyses  of  associated  with  a  lattice  F  and  an  ap¬ 

propriate  dilation  matrix  A.  In  Section  3  we  show  how  these  analyses  give  rise 
to  orthogonal  wavelet  bases.  Miscellaneous  remarks  and  acknowledgements 
are  collected  in  Section  4. 

We  use  standard  mathematical  terminology  and  notation.  A  brief  list  of 
some  of  the  conventions  used  here  follows:  The  Fourier  transform  /  of  an 
integrable  function  /  is  defined  by 

m  =  I  f{x)dx 

JR” 

and  distributionally  otherwise.  Basic  facts  concerning  Fourier  transforms  and 
distributions  will  be  used  without  further  elaboration  in  what  follows.  To 


avoid  the  pedantic  repetition  of  “almost  everywhere”  and  other  modifying 
phrases  which  axe  inevitably  necessary  when  dealing  with  functions  defined 
almost  everywhere,  all  equalities  between  functions  and  other  related  notions 
are  interpreted  in  the  distributional  sense  whenever  possible.  The  term  sup¬ 
port  is  also  used  in’ the  distributional  sense;  in  particular  the  support  of  a 
function  /  in  L\1R!^)  is  a  wdl  defined  closed  set.  If  W  is  a  collection  of 
tempered  distributions  then  WJs  the  collection  of  Fourier  transforms  of  el¬ 
ements  of  W,  in  other  words  W  =  {f  :  /  =  ^  for  some  g  in  W}.  For  a 
subset  Q  of  2R”,  a  linear  transformation  B  on  IR”,  and  an  element  y  of  2R" 
the  sets  B9.  axid  Q.  +  y  are  defined  by  BQ  =  {x  :  x  =  Boj  for  some  a;  in 
and  Q  +  y  =  {x  :  x  =  u  +  y  for  some  oj  in  fi};  L\^)  is  the  closure  of 
the  subspace  of  those  functions  in  L\JR!")  whose  support  is  contained  in  fi. 
Given  a  measurable  set  xn  denotes  its  characteristic  or  indicator  function 
and  |fi|  denotes  its  Lebesgue  measure. 


2  Multiresolution  Analyses 

2.1  Definitions 

Suppose  r  is  a  lattice  in  iR”,  that  is,  F  is  the  image  of  the  integer  lattice 
under  some  nonsingular  linear  transformation.  We  say  that  a  linear 
transformation  A  on  iR”  is  an  acceptable  dilation  for  F  if  it  satisfies  the 
following  properties: 

•  A  leaves  F  invariant.  In  other  words,  AF  C  F  where 
yir  ■=  {y  :  y  =  Ax  and  x  G  F}  . 

•  All  the  eigenvalues,  A,-,  of  A  satisfy  |A,|  >  1. 

These  properties  imply  that  |  det  A\  is  an  integer  a  which  is  >  2.  For  example, 
if  A  =  pi  where  p  is  an  integer  >  2  and  I  is  the  identity  then  A  is  an 
acceptable  dilation  for  any  lattice  F  and  a  =  |  det  A\  =  p". 

Such  an  A  induces  a  unitary  dilation  operator  Ua  •  f  — ^  UAf  on 
defined  by 

(3)  UaUx)  =  I  det  A\-^l'^f{A-^x)  . 
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If  V  is  a  subspace  of  L‘^{]R^)  we  use  the  customary  notation  UaV  to  denote 
the  image  of  V  under  Ua,  that  is,  UaV  =  {/  ■  f  =  UaQ,  9  €  V}.  The 
translation  operator  Ty  is  defined  by  Tyf(x)  =  f(x  -  y). 

A  multiresolution  analysis  V  associated  with  (F,  A)  is  a  family  {Vj}jQZ  of 
closed  subspaces  of  which  enjoys  the  following  properties: 

Al.  Vj  C  VjAi  for  all  j  in  Z. 

A2.  [Jj^^Vj  is  dense  in 

A3-  =  {0}- 

A4.  f{x)  e  Vj  if  and  only  if  f{Ax)  €  V^+i.  In  other  words 

Vj  =  f/^Vo,  j  e  2Z. 

A5.  Vo  is  invariant  under  r^.  More  specifically,  if  f{x)  is  in  Vo  then  so  is 
f{x  —  7)  for  all  7  in  r. 

A6.  There  is  a  function  0  €  Vb>  called  the  scaling  function,  such  that 
is  a  complete  orthonormal  basis  for  Vq. 

The  case  A  =  27  is  often  referred  to  as  a  dyadic  multiresolution  analysis. 
It  is  also  the  case  to  which  most  of  the  current  literature  is  devoted  and, 
except  for  certain  technicalities,  is  representative  of  the  general  case. 

A  multiscale  analysis  is  said  to  be  composed  of  generalized  spline  functions 
in  the  sense  of  Meyer  if  all  the  elements  of  the  subspace  V^o  are  continuous 
and  the  mapping  which  maps  /  into  the  sequence  of  values  {/(7)}7er  is  an 
isomorphism  from  Vq  onto  1^{T). 

Remark  The  above  definition  has  several  redundancies  and  thus  is  not 
very  compact,  for  example  see  [30].  However  we  use  it  here  for  reasons  of 
tradition  and  convenience. 

2.2  Lattices,  cosets,  and  related  items 

Every  lattice  T  in  JR”  is  the  image  of  ZZ"^  under  an  invertible  linear  trans¬ 
formation.  Thus  there  is  no  loss  of  generality  by  restricting  attention  to  the 
case  r  =  ZZ"^.  We  do  this  in  what  follows  to  avoid  unnecessary  obfuscation. 
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The  reader  who  is  interested  in  the  general  statements  of  the  results 
discussed  below  can  systematically  replace  2Z^  with  F  and  ‘I'kTZ^  with  the 
corresponding  dual  lattice  F'. 

Recall  that  if  A  is  an  acceptable  dilation  of  2Z^  then  ATZ'"  is  a  subgroup 
of  2^"  and  a  coset  of  ATL'"  is  a  set  of  the  form 

k  AA2Z^  =  {kA  Ay.3  e 

where  k  is  an  element  of  2Z^ .  Any  pair  of  cosets  are  either  identical  or  disjoint 
so  that  the  collection  of  all  cosets,  which  is  denoted  by  Z'^jAZ'^,  consists  of 
disjoint  sets  whose  union  is  2^".  The  number  of  disjoint  cosets  in  TZ'^fAZ'^ 
is  equal  to  a  =  |  det  A|.  A  subset  K.  of  is  said  to  be  a  full  collection  of 
representatives  of  if  it  contains  exactly  a  elements  and 

U  (/c  +  A^”)  = 

We  identify  the  representative  with  the  coset  it  represents.  For  instance, 
K  =  0  refers  to  the  element  in  A2^”  chosen  to  represent  this  coset. 

Since  the  Fourier  transform  of  f{A-^x)  is  af{A*^)  and  we  will  use  this 
fact  and  its  variants  often  in  what  follows,  we  adopt  the  notation  B  =  A*  to 
avoid  cumbersome  expressions  such  as  (A*)“^. 

Furthermore  we  write  )Ca  to  denote  a  full  collection  of  representatives  of 
2Z^/A2Z^, 

Since  B  is  an  acceptable  dilation  for  whenever  A  is,  we  write  ICb  to 
denote  a  full  collection  of  representatives  of 

■■■1  ^a}  ■ 

Observe  that  Ka  need  not  necessarily  be  a  full  collection  of  representatives 
of  and  vice  versa.  This  is  easily  seen  by  considering  the  two 

dimensional  example 

A  = 

Finally  we  note  that  any  absolutely  convergent  sum  over  can  be 
iterated  as  follows: 

(4)  E  =  E  I  E  • 

iez"  kg/Cb  Ogk+bz"  J 
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2.3  Scaling  functions 

An  application  of  Plancherel’s  formula  implies  that  the  orthonormality  of  the 
collection  -  j)}jez’'  is  equivalent  to 

(5)  E  =  1- 

kez" 

Thus  every  scaling  function  <f)  must  satisfy  (5). 

The  existence  of  a  scaling  function  is  equivalent  to  the  following  appar¬ 
ently  weaker  condition: 

A6i.  There  is  a  function  0  in  Vo  such  that  —  k)}kez’'  is  a  Riesz  basis 
for  Fo. 

In  many  specific  examples  this  property  is  much  more  apparent  and  easier 
to  verify  than  A6.  We  refer  to  the  function  (f)  of  condition  A6i  as  a.  pseudo- 
scaling  function. 

We  remind  the  reader  that  {4>(x  —  k)}kez’'  is  a  Riesz  basis  for  its  closed 
linear  span  if  and  only  if  there  are  positive  constants,  ci  and  C2  which  are 
independent  of  {ajtlfcez")  such  that 

Cl  E  /  „!  51  ak4>ix  -  k)\^dx  <  C2  E 

jfcez"  kez"  kez" 

It  is  well  known  and  easy  to  verify  that  {^{x  —  j)}jez’'  is  a  Riesz  basis 
for  the  closure  of  its  linear  span  if  and  only  if  there  are  positive  constants,  ci 
and  C2,  such  that 

(6)  Cl  <  E  <  C2. 

iez" 

Orthonormality  of  -  j)}jez’'  is  equivalent  to  the  case  ci  =  C2  =  1. 

Thus  the  alternate  item  A6i,  in  the  definition  of  multiresolution  analyses 
is  equivalent  to  saying  that  Vq  is  the  closed  linear  span  of  the  translates 
of  one  function  (f)  which  satisfies  (6).  The  original  item  A6  means  that  the 
constants  ci  and  C2  satisfy  ci  =  C2  =  1. 

2.3.1  Details 

The  calculational  technique  of  periodization  is  very  useful  in  the  study  of 
multiresolution  analyses.  It  can  be  described  as  follows: 
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Suppose  Q  is  a  tile  with  respect  to  a  lattice  T  for  iR".  That  is  Q  is  a 
compact  subset  of  iR”  which  enjoys  the  following  properties: 

•  U 

7er 

•  1(7  +  Q)  n  ^1  =  0  7  in  r  \  {0}. 

Then  if  h  is  any  intergrable  function  on  iR”  we  may  write 

(7)  j  h(i)di^jy.h(i+i)di. 

jjti 

In  the  applications  below  it  will  usually  suffice  to  take  Q  to  be  the  cube 
—  [-tTjTt]"  and  F  =  2'kZ'".  However  occassionally  other  variants  of  this 
formula  will  be  used.  Our  first  application  of  this  elementary  formula  will  be 
the  proof  of  (5). 

To  see  (5)  note  that  by  virtue  of  Plancherel’s  formula  ,  (7),  and  the  fact 
that  the  collection  is  orthonormal  we  may  write 

/  (l){x  -  j)<i){x  -  k)dx  =  f  =  Sj^k 

JR"  JQ^  m&Z" 

where  bj^k  is  the  Kronecker  delta.  The  last  equality  implies  (5). 

To  see  that  A6i  implies  A6  suppose  {<f>{x  —  k)}kez''  is  a  Riesz  basis  for 
Vq.  In  view  of  (6) 

_ 1 _ 

is  a  bounded  2'k2Z"'  periodic  function.  Since  /  is  in  Vb  if  and  only  if  it  enjoys 
the  representation 

m  =  mm 

for  some  F  in  L^(j!R"/27r.^”)  it  is  clear  that  the  function  0o  defined  by 

^  - 
y/Yljez”  ~  27i'j)P 

is  in  Vq.  Now  00  satisfies  (5)  so  the  collection  {(l>o{x  —  j)}jez’'  is  an  orthonor¬ 
mal  subset  of  Vq.  That  the  collection  {<l>o{x  -  j)}jez”  is  complete  follows 
from  the  fact  that  every  function  /  which  enjoys  the  above  representation  in 
terms  of  0  also  enjoys  an  analogous  representaion  in  terms  of  0o. 
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2.3.2  A  property  of  the  scaling  function  related  to  density 
Because  of  condition  A3  the  scaling  function  (j)  must  satisfy 

for  every  cube  Q  of  finite  diameter  in  A  detailed  explanation  of  this 
may  be  found,  for  example,  in  [30].  Thus  if  ^  is  continuous  at  the  origin  then 

(9)  |<^(0)1  =  1  . 

In  particular  if  (f>  is  integrable  over  iR”  then 

I  f  (j){x)dx\  =  1^(0)  1  =  1  . 

JR" 

2.4  Scaling  sequences 

In  view  of  A1  and  A4  the  scaling  function  <f>  must  satisfy  the  two  scale 
difference  equation 

(10)  <p{x)  =  Y,  Sk<t>{Ax  -  k) 

kez" 

for  some  sequence  {sjtlfcear"  The  sequence  {sfc}jtez"  is  called  the 

scaling  sequence.  The  Fourier  transform  of  (10)  is 

(11)  m  = 

where  B  =  A*  is  the  adjoint  of  A  and  S{^)  is  the  2'kZ”^  periodic  function 

5(0  =  -  E 

®  jfcez" 

2.4.1  Properties  of  scaling  sequences 

As  a  consequence  of  (5)  the  periodic  function  S  introduce  above  must  satisfy 

(12)  E  |S(«  -  2xB-'k)P  =  1 

Ke/Cfl 
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where  K.  is  any  full  collection  of  representatives  of  TZ^IBTZ,'^.  Furthermore 
if  is  continuous  at  the  origin  then  S{^)  is  also  continuous  at  the  origin 
and 

(13)  5(0)  =  1  . 

This  follows  from  (11)  and  the  fact  that,  in  view  of  (9),  (^(0)  7^  0. 

In  terms  of  the  scaling  sequence  (12)  is  equivalent  to 

(14)  S^—Aj^k  ~  0,^j,k  • 

kez" 

If  this  sequence  is  also  in  /^(^")  then  (13)  is  equivalent  to 

(15)  53  =  a  . 

kez'' 

Thus  (14)  and  (15)  are  a  necessary  set  of  conditions  for  an  sequence 

to  be  a  scaling  sequence. 

Details  To  see  (12)  use  (5),  (11),  iteration  (4),  the  fact  that  5(5"^^)  is 
2'kBZ'^  periodic,  and  (5)  again  to  write 

1  =  Ei.z-  \S{B-\^  -  2itmB-\i  -  2i^3)f 

=  E.€/Cb  {Ej€<+BZ-  -  27rK)  -  2irj)\^} 

=  -  27rK)P 

which  is  the  desired  result. 

2.4.2  More  properties  of  the  scaling  functions  and  sequences 

Note  that  in  view  of  (10)  0  may  be  considered  as  a  fixed  point  of  the  trans¬ 
formation  (f)  — >  Djtez"  Sk<p{Ax  —  k).  Unfortunately  the  solution  of  (10)  is  not 
unique  since  the  distribution  (j)Q  defined  by 

MO  =  mm 

is  also  a  solution  of  (10)  whenever  h  is  any  locally  integrable  function  which 
satisfies  h{B^)  =  h{i,).  Such  a  distribution  0o  may  fail  to  be  a  scaling  func¬ 
tion.  On  the  other  hand  if  h  also  enjoys  \h{^)\  =  1  for  almost  all  ^  then 
(j)o  is  also  a  scaling  function  but  not  necessarily  for  the  same  multiresolution 
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analysis.  Further  restrictions  are  needed  on  the  scaling  sequence  {sk}kez’' 
and  the  scaling  function  <f)  to  guarantee  a  unique  solution  of  equation  (10), 
for  example  see  [2,  13]  for  some  details  concerning  this  matter. 

The  scaling  function  (f)  is  not  unique  since  any  function  (f)i  whose  Fourier 
transform  satisfies 

MO = m)k() 

is  also  a  scaling  function  for  the  same  multiresolution  analysis  whenever  H  is 
any  measurable  2‘kZ’"  periodic  function  which  satisfies  \H{0\  =  1  for  almost 
all  C  The  scaling  equation  satisfied  by  <j>i  will,  in  general,  be  different  from 
the  one  satisfied  by  0. 

2.5  Examples 

2.5.1  Multiresolution  analyses  generated  by  self  similar  sets 

Suppose  K.A  is  a  full  collection  of  representatives  of  Z'^jAZZ'^.  Consider  the 
compact  set  Q  defined  by 

00 

(16)  Q  =  {xelR^:x  =  Y^  €,•  G  Ka)  • 

Note  that  Q  depends  both  on  A  and  the  choice  of  ‘digits’,  /C^.  This  set 
satisfies  many  intersting  properties.  For  our  purposes  it  suffices  to  note  the 
following: 

•  Q  is  self  similar  in  the  affine  sense.  That  is, 

(17)  AQ=  U  (^  +  ^) 

where  the  terms  in  the  union  are  essentially  mutually  disjoint  ,  that  is 
|(ki  +  Q)  n(«2  +  Q)1  =  0  whenever  ki  ^  K2- 

•  The  characteristic  function,  xq>  of  this  set  satisfies 

(18)  Xq(^)  =  E  XQ(^a:  -  k)  . 

k^Ka 
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This,  of  course,  is  equivalent  to 
{sfclfcez"  defined  by 


satisfies  properties  (14)  and  (15). 


(17).  Also  note  that  the  sequence 

if  A;  e  K,a 
otherwise 


•  The  measure  of  Q,  |Q|,  is  an  integer  >  1. 

The  characteristic  function  of  Q,  Xq>  generates  a  multiresolution  analysis 
in  the  following  sense:  If  Vq  is  the  L\SV)  closure  of  the  linear  span  of 
{Xg(a:  —  k)}kez”^  then  the  collection  of  spaces  V  =  {Vj}jez  defined  by 

=  Ua% 

=  {fix)  :  fiA-ix)  e  Vo} 


is  a  multiresolution  analysis  associated  to  i2Z^,A).  This  is  a  consequence  of 
the  fact  that  xg  satisfies  (18)  and  the  fact  that 

Z)  -  27rA:)|^  >  0  a.e. 

itez" 


See  [30]  for  more  details. 

Furthermore,  if  |Q|  =  1  then  xg  is  a  scaling  function  for  this  multireso¬ 
lution  analysis.  In  other  words,  {xg(^  “  k)}kez’'  is  a  complete  orthonormal 
system  for  Vq,  For  conditions  which  guarantee  that  |Q|  =  1  see  [17,  30]. 

On  the  other  hand  if  [Q]  >  1  then  {xg(^  —  k)}kez’'  fails  to  be  a  Riesz 
basis  for  Vq.  Nevertheless  it  is  not  difficult  to  find  a  scaling  function  for  V. 
The  function  <p  defined  by  the  formula  for  its  Fourier  transform 


m  = 


_ xc?(0 _ 

{Efcez"  |xq(^  ”  27rA:)p}^'' 


is  one  such  example.  In  many  cases  it  is  possible,  and  often  not  difficult,  to 
find  a  scaling  function  which  is  the  characteristic  function  of  an  appropriate 
set  with  measure  one;  however,  it  is  not  clear  whether  this  is  always  the  case. 

For  more  details  concerning  these  multiresolution  analyses  see  [17,  ?,  30]. 
We  conclude  this  subsection  with  several  specific  examples. 
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Univariate  examples  In  the  case  n  =  1  consider  A  =  2.  If  Ka  =  {0, 1} 
then  Q  is  the  interval  [0, 1]  and  xq  is  a  scaling  function  for  the  multiresolution 
analysis  which  it  generates;  note  that 

Vo  =  {f  £  L^{IR)  :  /  is  constant  on  the  intervals  {j,j  +  1),  j  €  Z)  . 

If  Ka  =  {0,  m}  where  m  is  an  odd  integer  ^  1  then  Q  is  the  interval  [0,  m] 
and  XQ  fails  to  be  a  scaling  function  for  the  multiresolution  analysis  which  it 
generates.  However  it  is  apparent  that  this  multiresolution  analysis,  namely 
the  one  generated  by  X[o,mi,  is  the  same  as  the  one  generated  by  X[o,i]- 

If  A  =  -2  and  tCA  =  {0, 1}  then  Q  is  the  interval  [-2/3, 1/3].  The  rest  of 
the  remarks  made  in  the  case  A  =  2  are  valid  mutatis  mutandis  in  this  case. 

The  case  A  =  3  is  more  interesting.  If  JCa  =  {0. 1)2}  then  Q  is  the  interval 
[0, 1]  and  Xq  is  a  scaling  function  for  the  multiresolution  analysis  which  it 
generates;  here  again  Vq  is  the  same  as  above.  If  1C  a  is  a  multiple  of  {0, 1, 2} 
then  Q  is  the  same  multiple  of  the  interval  [0, 1]  and  XQ  generates  the  same 
multiresolution  analysis  as  X[o,i]-  However  if,  for  instance,  ICa  =  {0,1,5} 
then  Q  is  a  disconnected  set  of  measure  one  and  the  multiresolution  analysis 
generated  by  Q  is  very  different  from  the  one  generated  by  [0, 1];  indeed  here 

Vo  =  {/  G  L‘^{IR)  :  f  is  constant  on  the  sets  j  +  Q,  j  E  ^}  . 

More  generally  one  may  consider  the  case  A  —  N  where  N  is  an  integer 
such  that  \N\  >  1.  The  cases  N  >  3  aie  not  unlike  the  case  N  =  3. 


Bivariate  examples  In  higher  dimensional  spaces  the  sets  Q  can  be  quite 
interesting.  We  will  limit  ourselves  to  several  standard  examples  in  IR^ 

First  consider  A  =  21  where  /  is  the  2  x  2  identity  matrix.  If  the  members 
of  ICa  a-re  described  by  the  columns  of 

/  0  1  0  1  \ 

\0  0  1  1  ) 


then  Q  is  simply  the  unit  square  [0, 1]^  and  xq  is  the  scaling  function  for  the 
multiresolution  analysis  which  it  generates.  However  if  the  members  of  ICa 
are  described  by  the  columns  of 


/  0  1  0  3  \ 
\  0  0  3  1  / 
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Figure  4:  The  set  Q  corresponding  to  A  =  21  and  JCa  described  by  the 
columns  of  (19) 


then  Q  is  set  described  by  the  shaded  area  in  Figure  4. 

If  / 

^  =  (-1  J 

and  the  members  of  K-a  are  described  by  the  columns  of 

0  1 
0  0 

then  Q  is  the  so-called  twin  dragon  of  appropriate  size  and  position  which  is 
described  by  the  shaded  area  in  Figure  5. 

For  more  examples  see  [17,  ?,  30]. 

2.5.2  Univariate  piecewise  polynomial  splines 

This  example  nicely  illustrates  the  usefulness  of  condition  A6i. 

Piecewise  linear  splines  First  consider  the  sequence  of  closed  subspaces 
of  L^(iR)  defined  as  follows: 

V)  =  {/  G  L^(iR)  :  /  is  continuous  on  2R 

and  linear  on  the  intervals  [k2~^,  {k  +  1)2“-'],  k  e  . 


Figure  5:  The  set  Q  corresponding  to  the  “twin  dragon” 

The  sequence  V  =  {Vj}j^z  is  a  family  of  closed  subspaces  of  L^{]R)  which 
satisfies  conditions  A1  through  A5  of  a  dyadic  multiresolution  analysis.  That 
condition  A6  is  satisfied  is  not  immediately  clear.  However  it  is  fairly  trans¬ 
parent  that  condition  A6i  is  satisfied.  Namely,  integer  translates  of  the 
function 
(20) 

constitute  a  Riesz  basis  for  ^0-  Thus  V  is  a  dyadic  multiresolution  analysis. 
The  scaling  function  whose  recipe  was  given  in  Subsection  2.3  is 

^-2 


and  the  scaling  sequence  is  the  sequence  of  coefficients  {sk}kez  of 


q(i:\  —  /  1^  +  47rfc| 

^^^~\Ekez\(  +  2^k\-^ 


-4'»l/2  1 


kez 


Note  that  because  of  the  analyticity  of  the  above  expressions  for  <j)o  and 
S{^)  both  4)o{x)  and  {sjfc}fcgz  enjoy  exponential  decay  as  x  and  k  go  to  ±oo 
respectively.  See  Figure  6. 
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N/ 


Figure  6:  The  piecewise  linear  pseudo-scaling  function  (f>  and  the  scaling 
function 

The  general  case  More  generally  we  may  consider  the  k  fold  convolution 
of  X,  the  characteristic  function  of  the  interval  [0, 1],  with  itself 

<f){x)  =  x*'y*xix)  • 
k 

This  function  may  also  be  described  via  its  Fourier  transform  as 


Note  that  (p  is  k  —  2  times  continuously  differentiable  and  coincides  with  a 
polynomial  of  degree  <  k  —  1  on  every  interval  {k,k  +  l),  k  E  2Z.  Also  observe 
that  in  the  case  fc  =  1  the  function  (p  is  simply  x  and  in  the  case  fc  =  2  the 
function  <p  is  same  as  the  function  described  by  (20)  translated  one  unit  to 
the  right. 

If  Vo  is  the  L^{1R)  closure  of  the  linear  span  of  {<p{x  —  m)}mez  then  the 
collection  of  spaces  V  =  {Vj}j^z  defined  by 

Vj  =  {fix)  :  /(2-ix)  €  Vo} 

is  a  dyadic  multiresolution  analysis  of  L^{JR).  As  in  the  case  of  the  mutires- 
olution  analyses  generated  by  tiles,  this  is  a  consequence  of  the  fact  that  <p 
satisfies 

•AW  =  2^  E  (  *  )  ■^(2^  -  “) 

"  m=0  \  / 
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and  the  faet  that 


a.e. 


E  l■^(^  -  2iri)P  >  0 

lez 

This  multiresolution  analysis  consists  of  piecewise  polynomial  splines  of  order 

k. 

It  should  be  mentioned  that  the  term  spline  is  used  here  for  traditional 
reasons.  In  the  case  that  k  is  an  even  integer  V  is  a  multiresolution  analysis 
composed  of  splines  in  the  sense  of  Meyer  whereas  in  the  case  that  k  is  an 
odd  integer  V  fails  to  be  such  a  multiresolution  analysis. 

The  subspace  Vq  is  the  intersection  of  L^(iR)  with  the  class  of  those 
tempered  distributions  s  whose  fc-th  order  derivative  is  a  distribution  of  order 
zero  supported  on  2Z,  in  other  words,  s  satisfies 


— 


(^)  =  £  -  m) 

mez 


where  6{x)  is  the  unit  Dirac  measure  at  the  origin.  This  description  of  Vq 
can,  of  course,  be  used  as  its  definition. 

As  in  the  piecewise  linear  case,  when  k  >  2  the  function  <f>  fails  to  be 
a  scaling  function.  On  the  other  hand  the  collection  {<f>{x  —  l)}i^z  is  a 
Riesz  basis  for  Vq,  and  the  function  (f>o  defined  by  the  formula  for  its  Fourier 
transform 

is  a  scaling  function.  More  details  concerning  this  multiresolution  analysis 
may  be  found  in  [5,  14,  33]. 


2.5.3  Multivariate  analogues  of  the  univeiriate  spline  examples 

There  are  many  generalizations  of  the  examples  in  the  previous  subsection. 
Three  classes  of  multivariate  analogues  which,  in  a  certain  sense,  reduce  to  the 
univariate  splines  considered  above  are  the  (i)  box  splines,  see  [22,  25,  34,  38], 
the  so-called  (ii)  polyharmonic  splines,  see  [28],  and  (iii)  the  multiresolution 
analyses  generated  by  the  k  fold  convolutional  iterates  of  the  characteristic 
functions  of  the  sets  described  in  Subsection  2.5.1,  see  [39].  Here  we  indicate 
a  few  details  concerning  the  last  two  classes. 
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Poly  harmonic  splines  We  begin  with  the  basic  setup. 

Let  SHk,  k  =  1,2,...,  be  the  class  of  those  tempered  distributions  s 
whose  fc-th  order  Laplacian  is  a  distribution  of  order  zero  supported  on 
in  other  words,  elements  s  in  SHk  satisfy 

A'‘s(x)  =  Ctn6{x  -  m) 

where  A  is  the  Laplacian  in  iR”  and  if  fc  >  2  then  =  AA*^“^. 

Next  we  recall  that  a  matrix  A  on  iR"  is  said  to  be  a  similarity  ifA  =  pAo 
where  Aq  is  an  orthogonal  matrix  and  p  is  a  real  number.  If  A  is  also  an 
acceptable  dilation  for  then,  of  course,  |det  A\  =  |p"|  is  an  integer  >  2. 
In  the  case  n  =  2  these  matrices  are  all  of  the  form 

or  f  ^ 

p  m  j  y  p  —m 

where  m  and  p  are  integers  such  that  +  p^  >  2. 

Finally  let 

Vo  =  SHkf]LHlR^) 

and  note  that  Vq  is  a  closed  subspace  of  L‘^{]R^)  which  is  not  {0}  only  when 
4A:  >  n.  Furthermore  if  4fc  >  n  the  collection  {^{x—my^rn^z” ,  where  function 
6  is  defined  by 

~  {Emez-  K  -  27rm|-4'=}i/2  ’ 
is  an  orthonormal  basis  for  Vq- 

Now,  if  A  is  both  a  similarity  and  an  acceptable  dilation  for  and 
B  =  A*  then 

(21)  =  53 

meZ” 

whose  Fourier  transform  is 


where 


5(0 


={ 


EmeZ"  1^  +  27rRm|  ^  -  5^ 


jm^Z 


n  +  27rm|“^*'  j  a 


meZ" 
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Note  that  because  of  the  analyticity  of  the  above  expressions  for  ^  and  S{C) 
both  ^(a;)  and  {sm}mez"  enjoy  exponential  decay  as  |x|  and  \m\  go  to  oo 
respectively. 

In  view  of  (21)  and  the  fact  that  {<f){x  -  is  an  orthonormal  basis 

for  Vq  it  follows  that,  if  A  is  both  a  similarity  and  an  acceptable  dilation  for 
then  the  collection  of  spaces  V  =  {Vj}jez  defined  by 

v;-  = 

=  Oi^)  ■  /(^  ^  ^0} 

is  a  multiresolution  analysis  associated  to  If  2k  >  n  then  this 

multiresolution  analysis  is  composed  of  splines  in  the  sense  of  Meyer.  For 
more  details  see  [28]. 

Note  that  in  the  case  n  =  1  these  classes  reduce  to  the  classes  of  piecewise 
polynomial  splines  considered  in  the  previous  subsection  which  are  of  even 
order. 

Multiresolution  analyses  generated  by  the  k  fold  convolutional  it¬ 
erates  of  certain  characteristic  functions  The  characteristic  functions 
are  of  sets  which  were  considered  in  Subsection  2.5.1.  The  basic  setup  is  the 
following: 

Suppose  fc  is  a  positive  integer,  A  is  an  acceptable  dilation  for  JCa 
is  a  full  collection  of  representatives  of  Z'^jAZ^,  and  Q  is  the  compact  set 
defined  by  (16).  Let  (f)  be  the  function  defined  by 

(j){x)  =  xi^) 

k 

where  y  is  the  characteristic  function  of  Q.  This  function  may  also  be  de¬ 
scribed  via  its  Fourier  transform  as 

m  =  mo}'  • 

Note  that 

(22)  (Pix)  =  Sm<t>{Ax  -  k) 

meZ” 

whose  Fourier  transform  is 
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where  B  =  A*  and 


s{i)  = 


1 

a 


Y.  . 

mez" 


The  function  (f){x)  has  compact  support  and  the  sequence  {sm}m6Z" 
a  finite  number  of  non- zero  terms. 

As  in  Subsection  2.5.1  the  function  (f>  generates  a  multiresolution  analysis 
in  the  following  sense;  If  Vo  is  the  L^{1R")  closure  of  the  linear  span  of 
{<f>{x  —  m)}mez’'  then  the  collection  of  spaces  V  =  {Vjjjez  defined  by 


1^-  = 

=  {/(^)  :  /(^  ^  ^o} 


is  a  multiresolution  analysis  associated  to  This  is  a  consequence  of 

the  fact  that  satisfies  (22)  and  the  fact  that 

—  27rm)p  >  0  a.e. 

m€Z" 


Furthermore,  if  |(5|  =  1  then  {4>{x  -  m)}mez”  is  a  Riesz  basis  for  Vq.  On 
the  other  hand  if  jQ]  >  1  then  {i){x  -  m)}mez''  fails  to  be  a  Riesz  basis  for 
Vq.  Nevertheless  in  either  case  it  is  not  difficult  to  find  a  scaling  function  for 
V.  The  function  0o  defined  by  the  formula  for  its  Fourier  transform 


m)  = 


m 

{Emez"  -  27rm)|2}^^^ 


does  the  job  in  both  cases.  Unfortunately  </>o  does  not  in  general  enjoy 
compact  support.  Other  interesting  facts  concerning  these  multiresolution 
analyses  can  be  found  in  [39] . 

Note  that  in  the  case  n  =  1  and  appropriate  choice  of  JC^  these  classes 
reduce  to  the  classes  of  piecewise  the  polynomial  splines  considered  in  the 
previous  subsection. 


2.5.4  Compactly  supported  scaling  functions 

Examples  of  compactly  supported  scaling  functions  can  be  found  in  Subsec¬ 
tion  2.5.1.  The  fact  that  the  spline  examples  considered  above  do  not  give 
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rise  to  compactly  supported  scaling  functions  suggests  that  examples  which 
are  smoother  than  those  considered  in  Subsection  2.5.1  may  be  more  difficult 
to  construct.  Indeed  the  construction  of  such  examples  in  the  univariate  case 
is  the  central  topic  of  Daubechies  work  [12],  see  also  [14]. 

The  simplest  example  considered  in  [12]  involves  the  function  (f)  described 
in  Example  1.2.3  and  Figure  2.  This  function  generates  a  dyadic  multiresolu¬ 
tion  analysis  of  L^{1R)  and  is  also  a  scaling  function  for  this  multiresolution 
analysis.  In  other  words,  if  Vq  is  the  L^(iR)  closure  of  the  linear  span  of 
{(f>(x  —  k)}kez  then  the  collection  of  spaces  V  =  {Vj}j^z  defined  by 

Vj  =  {fix)  :  fi2-^x)  e  Vo} 

is  a  dyadic  multiresolution  analysis  and  ^  is  a  scaling  function  for  this  mul- 
tiresolution  analysis. 

Recipes  for  constructing  dyadic  multiresolution  analyses  of  L^(iR)  with 
compactly  supported  scaling  functions  which  enjoy  any  desired  finite  order 
of  smoothness  can  be  found  in  [14]  along  with  many  specific  examples.  The 
rough  idea  is  to  construct  appropriate  scaling  sequences  and  define  the  cor¬ 
responding  scaling  functions  as  the  compactly  supported  solutions  of  the 
corresponding  two  scale  difference  equations.  The  difficult  part  is  guarantee¬ 
ing  that  the  resulting  solutions  are  indeed  scaling  functions  with  the  desired 
properties. 

Multivariate  examples  can  be  constructed  by  taking  appropriate  tensor 
products  of  univariate  examples.  For  example,  if  and  </>2  are  univariate 
scaling  functions  for  multiresolution  analyses  Vi  and  V2  which  are  associated 
to  the  dilations  x  — »  Nix  and  x  N2X  respectively,  may  or  may  not  be 
the  same  as  then  the  bivariate  function  (f)  defined  by 

(f)ix,y)  =  (t)iix)(j)2iy) 

is  a  scaling  function  for  a  multiresolution  analysis  of  L'^iJP?)  associated  to 
A)  where 

A  = 

On  the  other  hand,  in  the  case  of  a  general  acceptable  dilation  A  in  iR", 
n  >  2,  I  am  not  presently  aware  of  any  recipes  for  the  construction  of  com¬ 
pactly  supported  scaling  functions  for  multiresolution  analyses  associated  to 
i2Z'^,A)  which  enjoy  an  arbitrarily  high  finite  order  of  smoothness.  Some 
work  in  this  direction  may  be  found  in  [8]. 


24 


2.5.5  Band  limited  scaling  functions 

Suppose  n  is  a  compact  subset  of  iR”  which  enjoys  the  following  properties: 

•  C  BQ. 

•  |fl  +-27rk}\  =  0  for  any  element  fc  in  \  {0}. 

U  {Q  +  2Trk}  ~  . 

keZ’' 

•  For  all  cubes  Q  of  finite  diameter  in  IR" 

lim  r^TTT^  /  .  =  1  • 

j->oo  ^Q\  JB-iQ 

Consider  the  collection  of  subspaces  V  =  {Vj}jez  defined  by 
Vj  =  {fe  L\JBr-) :  supp/  C  . 

The  properties  of  Q.  imply  that  V  is  a  multiresolution  analysis  associated  to 
The  function  <p  defined  by 

(23)  kO  =  xkO 

is  a  scaling  function  for  this  multiresolution  analysis. 

Note  that  0  is  not  only  a  scaling  function  but  also  a  cardinal  function  for 
Vq.  Namely, 

0(fc)  =  5o,fc 

for  all  k  in  and  every  /  in  Vq  enjoys  the  representation 

/(^)  =  • 
kez” 

Since  the  last  equation  is  also  the  expansion  of  /  with  respect  to  the  orthog¬ 
onal  basis  {0(a;  -  fc)}jbez"  of  Vo  the  multiresolution  analysis  V  is  composed 
of  splines  in  the  sence  of  Meyer. 

Unfortunately  (j>{x)  has  poor  decay  properties  as  |a:|  goes  to  infinity.  In 
particular,  |a;|0(a:)  fails  to  be  in  L^(iR"). 

The  simplest  example  is  the  case  where  n  =  1,  is  the  interval  [— tt,  tt], 
and  A  =  B  =  2.  Here  0(a:)  =  sincx  and  we  refer  to  this  case  as  the  sine 
multiresolution  analysis. 
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Examples  with  better  decay  properties  Suppose  that  in  addition  to 
the  above  listed  properties  Q.  satisfies  the  following: 

•  The  interior  of  BQ  contains  Q. 

In  this  case  there  exit  functions  which  are  infinitely  differentiable,  are  iden¬ 
tically  equal  to  1  on  f),  and  have  support  in  Choose  one  which  is 

non-negative,  call  it  h,  and  define  (j)Q  via 

m = m  • 


Note  that  {4)q{x  -  k))kez^  is  a  Riesz  basis  for  the  closure  of  its  linear  span 
Vb.  If 

=  {f{x)  :  f(A-^x)  G  Vo} 

then  V  =  {Vj}jez  is  a  multiresolution  anlysis  associated  to  (^”,  A)  with  a 
scaling  function  (j)  defined  by 


HO 

{Lkez-  Hi  -  ■ 


Unlike  the  earlier  case,  <f)  is  not  a  cardinal  function  for  Vq.  Nevertheless  V  is 
composed  of  spline  functions  in  the  sense  oi  Meyer]  the  cardinal  function  A 
in  this  case  is  given  by 


HO 

12kez’'  H^  ~  27rA:) 


Furthermore  both  (j){x)  and  X{x)  decay  faster  than  the  reciprocal  of  any 
polynomial  as  |x|  goes  to  infinity;  in  other  word,  for  any  positive  integer  p 
both 

|x|^<?i)(a:)  and  |xl’’A(x) 

are  bounded  on  iR”. 

Unfortunately,  even  in  the  simplest  cases,  an  explicit  formula  for  either  (f> 
or  A  is  not  easily  available.  For  this  reason  in  certain  applications  one  may 
wish  to  sacrifice  some  of  the  decay  properties  by  choosing  less  regular  h. 
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3  Wavelets 

3.1  Introduction 

In  what  follows  A  is  an  acceptable  dilation  for  V  =  {Yj^jez  is  mul¬ 
tiresolution  analysis  of  associated  with  (^”,^4),  and  0  is  a  scaling 

function  for  V  =  {Vj}jez  with  corresponding  scaling  sequence  {sj}jez''-  Let 
b  =  a  —  1  where  a  =  |det  A|.  The  objective  of  Subsections  3. 1-3. 5  is  to 
outline  how  this  framework  gives  rise  to  a  wavelet  basis  and,  in  the  process, 
to  constructively  prove  the  following: 

Theorem  1  There  exist  b  functions  which  enjoy  the  following 

properties: 

•  If  ijjQ  =  (f)  is  a  scaling  function  for  V  then  the  collection 

{-lpj{x  -  k)}j^(ox...,b},k€Z’' 

is  a  complete  orthonormal  system  for  the  subspace  Vi . 

•  The  collection 


is  a  complete  orthonormal  system  for 

Remarks  The  functions  ,  •  •  • ,  ^6  are  a  collection  of  fundamental  wavelets 
associated  with  the  multiresolution  analysis  V.  The  development  below  will 
show  (i)  that  this  collection  is  not  unique  in  any  sense,  although  in  the  case 
0  =  2  there  is  in  some  sense  a  canonical  formula  for  obtaining  tfi  from  the 
scaling  function  ifo,  and  (ii)  that  b  is  the  minimal  number  in  this  collection. 

An  important  objective  is  the  construction  of  wavelets  ipi,. . .  .ipb  which 
have  the  same  decay,  for  large  [a;],  that  the  “best”  scaling  function  has.  For 
example  if  V  has  a  compactly  supported  scaling  function  then  one  would 
like  the  if's  to  have  compact  support.  The  construction  outlined  below  does 
not  in  general  guarantee  this.  However,  because  of  the  importance  of  this 
matter,  we  will  make  remarks  concerning  this  issue  at  convenient  places  in 
the  development. 
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3.2  The  basic  setup 

Let  Wj  be  the  orthogonal  complement  of  Vj  in  In  other  words 
(24)  Vj+i  =  Wj  ©  V;- 

where  ©  denotes  the  orthogonal  direct  sum  of  linear  subspaces.  Iterating 

(24)  results  in 

Vj+i  =  Wj  ©  Wj-i  ©  •  •  •  ©  Wj-k  ©  Vj-k 

for  any  positive  integer  k  The  last  identity  together  with  properties  A2  and 
A3  imply  that 

(25)  L\1R^)  =  E  ®Wj 

jez 

Now,  the  fact  that  mapping  Ua  defined  by  UAf{x)  =  a~^/^/(A~^x)  is  a 
unitary  transformation  together  with  property  A5  imply  that 

(26)  Wj  =  UJ^'Wo 

for  all  integers  j.  Thus  knowledge  of  the  structure  of  the  subspace  Wq  gives 
us  analogous  knowledge  of  Wj  for  all  j.  In  particular,  if  there  are  functions 
^1, ...,  V’6  such  that 

(27)  -  fc)},g{i,...,6},fcez" 
is  a  complete  orthonormal  system  for  Wq  then 

is  a  complete  orthonormal  system  for  Wj. 

In  what  follows  we  will  constructively  show  that  there  are  6  =  a  —  1 
functions  V’l,  •  •  • ,  ^”6  whose  Z"  translates  (27)  are  a  complete  orthonormal 
system  for  Wq.  In  view  of  (25)  and  (26)  we  may  conclude  that 

(28) 

is  a  wavelet  basis  for  L‘^{1R^). 

Since  Wq  C  Vi  it  follows  that  the  fundamental  wavelets  •0i,  •  •  ■ ,  ^>6  satisfy 
(29)  ^pi{x)  =  Y,  Si,k(f>{Ax  -  k) 

kez" 
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for  i  €  {1, . . . ,  6}  where  0  is  a  scaling  function  for  the  multiresolution  analysis 
V  and  are  appropriate  sequences  in  Now,  if  {so,fc}fceZ"  = 

{sfc}jtGZ"  is  the  scaling  sequence  for  <f),  namely 

4>{x)  =  51  ^Q,k(f>{Ax  -  k)  , 
kez” 

then  from  the  orthogonality  relations 

I  (x  -  fci ) {x-k2)dx  =  ,  ,.-2  6k^  ,fc2 

JR" 

where  V’o  =  *i  and  Z2  are  indicies  taking  values  in  and  ki 

and  k2  are  elements  in  we  may  conclude  that  the  sequences 
i  =  0, 1, . . . ,  m,  must  satisfy 

(30)  53  ^h,k-Aj^h,k  ~ 

kez" 

for  j  G  2Z"^. 

Conversely,  if  {si^k}k€Z",  i  =  are  sequences  in  P{2Z’^)  which 

satisfy  (30)  where  if  {5o,fc}fceZ"  is  the  scaling  sequence  for  (j)  then  the  functions 
^>1, . . . ,  V’m  defined  via  (29)  are  candidates  for  fundamental  wavelets.  Indeed 
the  constructive  proof  of  the  existence  of  a  collection  of  fundamental  wavelets 
essentially  relies  on  the  construction  of  such  sequences.  The  details  however 
are  carried  out  in  the  Fourier  transform  or  frequency  domain. 

3.3  A  characterization  of  Wq 

Recall  that  ^  is  a  scaling  function  for  V  =  with  corresponding  scaling 

sequence  {sj}j^z"  a^nd 

m  =  s{B-'mB-'o 

where  B  =  A*  is  the  adjoint  of  A  and  S{^)  is  the  27rZ”  periodic  function 

S(0  =  i  j: 

®  keZ" 

Also  recall  that  )Ca  and  K-b  denote  full  collections  of  representatives  of 
JAZ'^  and  Z”^  jBZ'^  respectively.  See  Subsection  2.4. 

Note  that  every  2'kBZ’^  periodic  function  is  of  the  form  F{B~^^)  where 
F{^)  is  2'kZ'^  periodic. 
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Proposition  1  The  function  f  is  in  Wq  if  and  only  if  f  enjoys  the  repre¬ 
sentation 

(31)  /({)  =  F{B-'OkB-'0 

where  F{^)  is  in  12'kZ'")  and  satisfies 

(32)  E  -  2nu))S{B--^{C-2'Ku))  =  0 

velCa 

for  almost  all 

Proof:  Observe  that  every  element  p  in  Vq  enjoys  the  representation 

(33)  m  =  Giom 

where  G{^)  is  in  L^(JR"/27r^")  and,  conversely,  every  such  G  gives  rise  to 
an  element  5'  in  Vb  via  (33).  Since  gl^)  is  in  Vo  if  oiiiy  if  9iB  i®  i^ 
Vi  and  Wo  C  Vi  we  may  conclude  that  reprsentation  (31)  follows  from  (33). 

Relation  (32)  is  equivalent  to  fact  that  /  is  orthogonal  to  Vq.  To  see  this 
use  the  following: 

•  Plancherel’s  formula, 

•  the  Fourier  transform  of  the  scaling  relation  (11), 

•  the  periodization  trick  (9),  and 

•  the  fact  that  if  H{^)  is  a  27r.^”  periodic  function  then 

Y,  =  E  -  2n„))  . 

jez"  v£Kb 

This  last  item  follows  from  (5),  the  periodicity  of  H,  and  an  application 
of  (4). 

Explicitly,  for  every  g  in  Vq  write 

f  f{x)g{x)dx= 

J IV^  J  R 

=  /  F(B-H)S(B-H)G{0\kB-^OM 

JR" 
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(vEICb  ) 


or,  more  succinctly, 


(34)  [  f{x)g{x)dx  = 

J 


/  Y,  F(B-'(f-27ri/))S(B-'(f-2jri/))  G(B{K. 

{i/eiCB  ) 

Thus,  if  (32)  holds  then  (34)  implies  that  /  is  orthogonal  to  all  g  in  Vo- 
Conversely,  if  /  is  orthogonal  to  all  ^  in  Vo  the  integral  on  the  right  hand 
side  of  (34)  must  vanish  for  all  G  in  L\JR!'I2'kZ,^)  which  implies  (32).  ■ 


3.3.1  An  alternate  expression  for  T,F{B  ^(^  -  2'kv))S[B  ^(^  -  2'kp) 
The  27r.^"  periodization,  or  average  over  the  group  of 

F{B-^i)S{B-^0 

on  the  left  hand  side  of  formula  (32)  suggests  a  kind  of  inner  product  which 
maps  pairs  of  2TrB2Z^  periodic  functions  into  a  2TrZ'^  periodic  function. 
This  “product”  may  be  more  conveniently  expressed  in  terms  of  the  27: 
periodic  components  of  the  respective  27: BZ^  periodic  functions.  To  wit: 
Every  27:BZ'"  periodic  function  F{B~^^)  enjoys  the  representation 

(36)  F(B-H)  =  4=  E 

V  ®  k^ICa 

where  Fk(^),  /c  G  is  a  27: Z"^  periodic  function  which  can  be  derived  from 

F{B~^^)  via  the  formula 

(36)  B.(0  =  4=  E  ns-'n  -  27r^))e<'‘-‘“’«-^'-l  . 

V  ®  j/€/Cb 

If  G{B~^C)  is  another  27: BZ""  periodic  function  then 

(37)  Y.  Y  FM)GM 

I/£/Cb  KEICa 
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where,  of  course,  is  related  to  G  in  the  same  way  that  Fk  related  to  F. 
Note  that  (37)  gives  us  an  alternate  expression  for  the  inner  product  defined 
by  the  right  hand  side  of  (32). 

Formulas  (32),  (35),  and  (37)  suggest  that  2TrBZ^  periodic  functions 
F{B~^0  be  viewed  as  a  tuples  of  27r2^”  periodic  functions 

where  (/ci, . . . ,  Ko)  is  an  ordering  of  ICa  and  , . . . ,  F^^  are  related  to  F{B~^^) 
via  (36).  We  do  this  in  what  follows. 


3.3.2  Details 

We  begin  by  noting  the  following: 

Lemma  1  If  u  £  Kb  then 

(38)  53  g.-2tr<A->K,./>  _  _ 


This  is  nothing  more  than  the  observation  that  the  mapping 


K  + 


is  a  character  of  the  coset  group  and  (38)  is  a  well  known  relation 

for  such  functions.  For  example  see  [36].  For  the  sake  of  completeness  we 
recall  the  argument. 

Proof:  If  1/  =  0  then  (38)  is  transparent,  li  u  ^  0  then  there  is  a  yii  in 

Ka  such  that  ^  l.  Since  Ka-  fiis  also  a  full  collection  of  coset 

representatives  so  that 


k&Ka 


_  gt27r(A-*K,i/) 


? 


we  may  write 


pi2'rr(A  ^k,u) 


_  g«2x(A-‘(#i,i/)  Qi2ir{A  ^{k-b),i') 

_  g<2x(A-V.''>  ^2it{A-^K,v) 


Since  #  1  the  last  string  of  equalities  implies  (38). 
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Formulas  (35)  and  (37)  are  easy  consequences  of  this  lemma.  To  see  (35) 
use  identity  (36)  for  Fk,  interchange  order  of  summation,  and  apply  formula 
(38)  to  write 


4=  E 

V  ®  kg/Ca 


«g/Ca  i  ^  I'eK.B  J 

=  -  F(5-^(^-27ri/))|  53  I  =  F(B-^0  • 

^  u&Kb  (<cGACa  J 

To  see  (37)  use  formula  (36)  for  and  interchange  various  orders  of 
summation  appropriately,  and  use  identity  (38)  to  write 


E  m)GM 

kG/Ca 

kG/Ca  i.  i/G/Cb  j 

{-4  E  G'(B-He-27rAx))V^"''’^-2"'^4 

J 


=  E  E  n5-He-27ri/))G(B-He-27r/i))x 

v^K-b 


j  1  v  --•27r(A-i/c,./-Ai)  I 

V^A  J 


=  E  F(B-'(e-27ri/))G(B-He-27ri/)) 

I'G/Cb 


Remark  We  alert  the  reader  familiar  with  the  application  of  these  formulas 
in  [25]  that  a  slightly  different  normalization,  which  results  in  a  factor  of  a 
in  the  right  hand  side  of  formula  (37),  was  used  there. 
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3.4  Consequences  of  the  characterization  of  Wq 

In  what  follows  we  use  the  following  conventions: 

•  0  =  1  det  A\  and  b  =  a  —  1 

•  {ko,  . .  ■ ,  and  {i/q,  .  ■  • ,  are  fixed  orderings  of  JCa  and  ICb  respec¬ 
tively  with  Ko  =  1^0  =  0- 

•  If  is  a  27r^"  periodic  function  then  F{C)  =  {FoiO^  -  ■  ->  ^biO) 

where 

V  ®  ueKa 

for  I  e  {0,...,b}.  We  identify  with  F{B-'^0  ^  indicated  in  the 
previous  subsection,  namely, 

(39)  =  • 

V  ^  /=0 

•  If  F(B~^^)  and  are  2'kB2Z’"  periodic  functions  then 

(F(0M))B='tFi{()W)  • 

/=0 

•  The  class  of  all  measurable  27r^"  periodic  vector  fields 

^(0  =  we),..-,n(0) 

with  values  in  (27“  is  denoted  by  yF(i?"/27r^”,(r“).  The  subclass  con¬ 
sisting  of  those  elements  of  yF(F"/27r2^",  (r“)  which  are  locally  square 
integrable  is  denoted  by  L\VF{R^  12'kZZ'^ ,(C^))  and  for  such  elements 

(F,G)=  /  {F{0,G{0)Bd^  . 

*'Qir 

In  view  of  the  observations  made  in  Subsection  3.3.1  the  conventions 
allow  us  to  express  the  identity  (12)  involving  the  scaling  factor  S{B 
and  identity  (32)  characterizing  the  periodic  factors  in  the  Fourier  transform 


34 


of  Wo  in  a  more  convenient  way.  Namely,  if  is  a  scaling  factor  then 

the  corresponding  element  S(^)  in  must  satisfy 

(40)  (5(0,5(0)b  =  1 

and  if  is  the  Fourier  transform  of  /  in  Wq  then 

(41)  {F{OJiO)B  =  G 

for  almost  all  Note  that  identity  (40)  implies  that  5(0  is  a  27r^"  periodic 
vector  field  with  values  on  the  unit  sphere  in  (F".  Relation  (41)  means  that 
f{()  is  a  27r^"  periodic  vector  field  which  is  orthogonal  to  5(0- 

Recall  that  the  Fourier  transform  of  every  element  in  Vi  may  be  expressed 
as 

/(o  =  F{B-^oko 

where  F(0  is  in  I2'k2Z'^).  In  view  of  the  identification  of  F(J5"^0 

with  F(0  we  see  that  every  /  in  Vi  can  be  identified  with  a  unique  F(0  in 
L2(yF’(i?"/27r^”,<r“))  in  a  natural  way. 

Suppose  Xo,  •  •  •  is  a  basis  for  FF(R"/27r.^",(r'“)  in  the  sense  that 
the  Gramm  matrix 

(42)  Gr(Xo,  .  .  .  jXft)  =  Xm{^)) B)l,m=0,...,b 

is  both  uniformly  bounded  and  uniformly  strictly  positive  definite,  that  is, 
there  are  positive  constants  c  and  C  such  that 

(43)  c  Yu  |r,„P  <  Z]  13  <CY 

m=0  /=0  m=0  m=0 

for  almost  all  ^  and  all  complex  a  tuples  (tq,  . . .  Then  every  F(^)  in 
L^(VF(jR”/27r^",(F“))  can  be  expressed  as 

1=0 


for  some  uniquely  determined  element  B(^)  =  (Bo(^), B'b(O)  in 

L\VF(R^/2%ZS^,(D^)). 
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Consider  the  functions  il^oix), . . . ,  defined  by 

(44)  *({)  =  X,(B-H)kB-'e) 

where  Xq,  . . .  ,^6  are  a  basis  of  the  type  discussed  above.  We  remind  the 
reader  that  </>  is  the  scaling  function  and 

=  4=  E 

V  ^  m=0 

where  Xi{0  =  •  •  ■  >  I  e{0,...,b}-  In  view  of  the  observa¬ 

tions  made  in  the  previous  paragraph  the  2Z'^  translates  of  the  ^/’s,  namely 

constitute  a  Riesz  basis  for  Vi . 

In  paricular  if 

5'o(0)b  =  0  for  all  ^  and  I  £  b}, 

then 

is  a  Riesz  basis  for  Wq.  In  analogy  with  the  notion  of  pseudo-scaling  func¬ 
tion,  we  refer  to  the  members  of  such  a  basis  as  pseudo-wavelets  and  the  set 
•0i(a;), . . . ,  as  a  full  collection  of  fundamental  pseudo-wavelets. 

Furthermore  if  Xq  =  5  and 

(X/(^), Xk{0)B  —  bi,k  for  3.11  ^  and  l,k  £  {0,  ■  ■  ■  ,b}  , 

then 

{i^i{x  -  j)}ie{i,...,b},jez" 

is  an  orthogonal  basis  for  Wq.  In  other  words,  tl;i{x), . . .  is  a  full 

collection  of  fundamental  wavelets. 

We  summarize  some  of  these  observations  as  follows: 

Proposition  2  Suppose  Xq,  . . .  ,X{,  is  a  basis  for  V^P(i2"/27r^  ,(C  )  in  the 
sense  that  the  Gramm  matrix  (42)  satisfies  (43).  Then  the  translates  of 
the  functions  iPq{x),  . . . ,  ipbi^)  defined  by  (44) j  namely  the  functions  of  x 

-  j)}ie{o,...,b},jez”  , 
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are  a  Riesz  basis  for  Vy .  In  particular  if  Xq  =  S  and 

(Xi(0,  ^fc(O) B  =  deltai^k  for  all  ^  and  A:  G  {0, . . . ,  6}  , 


then 

-  j)}ie{h...,b},jez'' 

is  an  orthogonal  basis  for  Wo  ond  thus  ^’1(^)5  •  •  • » ®  f'O'll  collection  of 
fundamental  wavelets. 


3.5  A  recipe 

It  is  now  clear,  in  principle  at  least,  how  to  obtain  a  full  collection  of  fun¬ 
damental  wavelets  in  terms  of  a  scaling  function  </>(a:)  and  the  corresponding 
scaling  factor  which  are  associated  with  the  multiresolution  analysis 

V.  The  recipe  goes  as  follows: 

•  First,  for  notational  convenience,  set  =  S{B~^^).  Next,  iden¬ 

tify  So{B~^^)  with  the  vector  field 


5o(e)  =  («),•••,  MO) 

in  (i2’"/27r^”,(r“).  The  components  of  5o(0  can  be  evaluated  in 
terms  of  Sq{B~^^)  via  formula  (39) 

•  Select  b  elements 


S„{()  =  (5™o«).  ■  •  ■ .  S„i,m  ,  m  e  {1. . 

in  Vf  (B"/2?rZ",lC‘‘)  such  that 


(45) 


/"SooCft  SoiiO  S»{e)\ 
Sio(«  SniO  ■■■  Su(0 

I  SboiO  ■■■  Sth(^)  j 


is  a  unitary  matrix  for  all  0  la  other  words,  Si,...,Sb  should  be  chosen 
so  that 


{Si{0,Smi0)B  =  S,, 
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for  all  ^  and  for  all  I  and  m  in  {0, ,  b}.  This  should  be  possible,  in 
principle,  since 

{So{OJoiO)B  =  l 

for  all 

•  Define  ^/,  /  =  1, . . . ,  6,  via  the  formula  for  its  Fourier  transform 

(46)  Uo  = 

where  ^  is  the  scaling  function  associated  with  the  scaling  factor  So{B~^^), 

m=0 

and  {/Co, . . . ,  K6}  is  a  full  collectiion  of  representatives  of  the  coset  group 
Equivalently 

=  E  -  j) 

j&Z" 

where  the  sequence  {sij}j^z”  is  the  sequence  of  coeflBcients  in  the  ex¬ 
pansion 

Si{()  =  i  E 

“  iez" 

By  virtue  of  Proposition  2  the  collection  , . . . ,  V’6  defined  by  (46)  is  a 
full  collection  of  fundamental  wavelets. 

The  only  unclear  instruction  in  this  procedure  involves  the  selection  of 
elements  5i, . . . ,  ^6  in  yF’(iZ”/27r.^”,(r“)  so  that  (45)  is  a  unitary  matrix. 
We  clarify  this  by  first  considering  specific  cases. 

3.5.1  The  case  a  =  2 

In  this  case  we  don’t  have  much  choice.  The  selection 

5i(0  =  (5^,-5^) 

does  the  job  nicely.  Indeed,  every  possible  choice  is  of  the  form 

5i(0  =  n0(5^,-5^) 

where  T(^)  is  a  measurable  27r.^”  periodic  function  which  satisfies  \T{^)\  =  1 
for  all 
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3.5.2  The  case  a  =  3 

Once  Si  is  chosen  then  S2  must  be  of  the  form 

52(0  =  r(0(5o(0  X  5i(0) 

=  T(^^)(^SqiSi2  —  5'o2'S^11)  ^02^10  “  SqqSi2^  SqqSii  ^Ol^io) 

where  T(0  is  a  measurable  27r^”  periodic  function  which  satisfies  \T(^)\  =  1 
for  all  ^  and  x  denotes  the  cross  product  in  which  is  essentially  defined 
by  the  second  equality. 

This  is  an  easy  consequence  of  the  fact  that  any  three  vectors  u,  v,  and 
w  in  satisfy  the  identities 


f  w\ 


{w,  u  X  v)  =  det 


u-j 


and  X  • 


The  first  implies  that  S2  is  orthogonal  to  both  So  and  Si  while  the  second 
implies  that  (52(05  52(0)b  =  1 

The  choice  of  5i(0  is  not  so  clear  however.  On  the  other  hand  if  the  first 
component  of  So{0,  namely  the  function  5oo(0>  enjoys  the  property  that  for 
some  positive  constant  c 

(47)  l-5oo(01  >  e 

for  all  ^  then  the  choice 

g  ^  (3^,-^,0)  (3^,0,-^,) 

\/l5oop  +  |5oip  \/|5oop  +  |5o2p 

will  do  the  job.  Analogous  choices  should  be  clear  if  5oi(0  or  5o2(0  satisfy 
(47).  If  none  of  the  components  of  5o(C)  is  uniformly  bounded  away  from  0 
then  an  appropriate  combination  will  to  the  job.  For  instance, 

Si  =  Ti(^)(5oi,  — 5oo,0) +  T2(^)(5o2,0,  — 5oo)  • 

Note  that  Ti{^)  and  T2(0  can  be  chosen  so  that  5i(0  is  orthogonal  to  5o(0, 
satisfies  the  identity  (5i(^),  Si{^))b  =  1  for  all  and  is  almost  as  smooth  as 

5o(0- 
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3.5.3  The  general  case:  part  1 


For  simplicity  write 

(48)  So(0  =  iMO,---MO) 

and  assume  that  one  of  its  components,  which  without  loss  of  generality  we 
take  to  be  UQ)  satisfies 


(49)  |«o(0l  >  c  >  0 

for  all  Define  Xi  =  (Xio, . . . ,  =  (Xm,  . Xbt)  by 


(50) 


f  Xr  \ 

/  —ui  uq  0  •  •  •  O'' 

X2 

—  «2  0  Mo  •  •  •  0 

♦  •  .  ♦  •  ♦ 

KXb) 

\  — M5  0  0  •••  Mo  ; 

and  observe  that 

(51)  (X,(OJo(0)b  =  0 

for  all  ^  and  all  I  in  For  each  I  in  the  vector  field 

Xi  =  (Xio, Xib)  may  also  be  described  by 


f  -M/ 


X,^  =  { 


Uo 

0 


if  m  =  0 
if  m  =  / 
otherwise. 


Next  note  that  the  Gram  matrix 

Gr(Xi,  ...,X,)  =  {{XliO,XmiO)B)l,m=l,...,b 


can  be  expressed  as 

(52)  Gr(Xi,...,X6)  =  |«o|''l6  +  H"P 


where  is  the  6x6  identity  matrix, 

b 

u  =  {ui,...,ub)  and  = 

j=i 
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The  matrix  P  is  the  orthogonal  projection  onto  the  one  dimensional  subspace 
generated  by 


so  that  in  matrix  notation 

|M|2p  =  u*u  . 

Relations  (51)  and  (52)  imply  that  the  Gram  matrix 

Gr(/So,  ^1,  — , 

is  uniformly  strictly  positive  definite,  in  other  words  it  satisfies  (43),  and 
that  the  functions  , . . . ,  V’t  defined  by  (44)  are  a  full  collection  of  pseudo 

wavelets.  ■ 

Orthonormal  vector  fields  Si, ...  ySh  which  are  orthogonal  to  5'o  can  now 
be  obtained  by  applying  some  sort  of  orthonormalization  procedure  to  the  col¬ 
lection  Xi,...,Xb.  For  example  one  may_apply  the  cl^sical  Gram  Schmidt 
method  to  some  permutation  of  Xi, . . . ,  to  obtain  S\,... ,  Si- 

A  more  symmetric  method  of  obtaining  Si, . . .  ,Sb  goes  as  follows:  First 

write  ,  _ 

Xi 

X2 

G  =  Gr(Xi,  ...,Xb)  and  X  =  . 

Now  recall  that  G  =  XX*.  Since  G  is  uniformly  strictly  positive  definite  we 
may  compute  both  G“^  and  its  square  root,  G  Let 

(53)  S  =  G-'/^X 

where 

S  = 


41 


Observe  that 


SS*  =  =  G-^/^GG“^/2  ^ 

which  implies  that  the  rows  of  S  are  orthonormal. 

The  method  for  obtaining  5i, . . . ,  5^  which  is  succinctly  summarized  by 

(53)  gives  rise  to  explicit  expressions  for  these  vector  fields  because  G  can 
be  expressed  as 

(54)  G  =  (|«oP  +  1«P)P  +  l«ol^(l6  -  P)  , 

which  is  simply  a  convenient  way  of  writing  equation  (52).  From  (54)  it 
follows  that 

G-'/2  ^  ^  |„[2)-l/2p  +  _  p)  , 

which  may  be  re-written  as 

= ((i«op + -  i“or‘)  p + Wor'ii 

or,  which  is  even  more  convenient,  as 


f  -«i 

j3uiui  -1-  a 

/3uiU2 

G-^/^X  = 

-U2 

pU2U2  -1-  a 

; 

where 

\  -Ub 

PUbUi 

a  = 

PUbU2 

«o 

hol 

pUiUb 

f3u2Ub 

(3ubUb  -f-  a 


and 


/?  = 


—a 

1  +  l«o|  ' 
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In  other  words,  in  this  case  for  each  I  in  6}  the  vector  field  Si  — 

(5/0, . . . ,  Sib)  may  be  described  by 


(56) 


Slm=< 


-Ul 

f3uiui  +  a 
pUiUm 


if  m  =  0 
if  m  =  Z 
otherwise. 


We  summarize  this  discussion  as  follows: 

Proposition  3  If  one  component  of  Sq  is  uniformly  bounded  away  from  0 
then  the  above  procedures  give  rise  to  explicit  formulas  for  mutually  orthonor¬ 
mal  vector  fields  5i, . . . ,  56  in  I/P(il”/27r.S”,(C'“)  which  are  orthogonal  to  Sq. 
The  corresponding  functions  V’l,  —  ,tpb  defined  by  (44)  a  full  collection 
of  orthonormal  wavelets. 

In  particular,  if  So(^)  =  . . .  ,Ub{C))  and  lno(0l  ^  ^  ^  ,  ^ 

the  vector  fields  Si,...,Sb  whose  components  are  described  by  (56)  are  mu¬ 
tually  orthonormal  and  are  orthogonal  to  Sq.  The  corresponding  functions 
ipi,...,‘ipb  defined  by  (44)  are  a  full  collection  of  orthonormal  wavelets. 


Remark  There  are  many  interesting  cases  where  the  hypothesis  of  the 
above  proposition  is  valid.  For  instance,  suppose  the  scaling  factor  S{B 
is  real  valued  and  non-negative;  examples  of  this  will  be  given  in  Subsection 
3.7.  The  first  component  uq  of  5o(^)  is  given  by 

“0  =  4=  E 


Since 

1=  E  |5(B-‘(f  -  2»-./))p  <  E  5(S-'K-2x./))<a 

i'G/Cb 

it  follows  that  ^ 

-p  <  no  <  \/a  • 

V« 

Details  Note  that  in  matrix  notation  X  can  be  expressed  as 

X  =  [«*,nol6] 
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and  hence 


G  =  XX*  =  K,ttol6] 


—u 

Uolb 


=  U*U  +  llioplfc  • 


Thus 


g-1/2  ^ 


Yl«oP  +  N"  l“ol/  I'' 

and,  since  |uoP  +  |up  =  1,  we  may  write 


1  \  1  »  1  T 

-  1 — r  I  r-rT«  u  + - It 

2  uo 


G-'/^x  = 


'  1  \  1  ,  1  ^ 

^  |«o|  j  l^^ol 


-?X*,«ol6] 


-U 


■HCiir'"") 


Using  the  substitution  [up  =  1  —  |noP  =  (1  ~  ko|)(l  +  l^ol)  in  tiie  Inst 
expression  results  in  (56). 


3.5.4  The  general  case:  part  2 

Finally  suppose  no  components  of  So{^)  =  (ui(^),...  ,Ub{^))  are  uniformly 
bounded  away  from  zero.  Since  {So,Sq)b  =  1  it  follows  that  there  are  mea¬ 
surable  sets  Qo,...,^b  which  enjoy  the  following  properties 

•  Each  set  is  27r.^"  periodic.  In  other  words,  if  ^  is  in  Vim  then  ^  -1-  27rfc 
is  in  for  all  k  in  and  this  hold  for  all  m  in  {0, ... ,  b}. 

•  Whenever  ^  is  in  Q,rn  then 

(57)  WmiOl  >  ^  • 

•  They  are  mutually  disjoint.  That  is,  flj  fl  =  0- 

•  Their  union  covers  iR".  That  is, 

(58)  U  =  iR”  . 

Tne{0,...,6} 
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Let  Xo(0»  •  •  •  >  XbiO  Le  the  characteristic  functions  of  of  fio, .  •  • ,  respec¬ 
tively.  The  last  two  listed  properties  of  the  fl’s  imply  that 

^^(0  =  1 

m=0 

for  all  In  view  of  the  property  described  by  equation  (57)  using  the  meth¬ 
ods  outlined  above  for  each  m  in  {0, . . . ,  6}  we  can  construct  vector  fields 
,  5™  in  V’F(i2”/27r.^”,(r“)  such  that 

{S]^Jo)b  =  0  for  all  fc  6  6} 

and 

(5r,  S^B  =  Sk,iXm{0  for  all  fc,  1  e  {1, . . . ,  6}  . 

These  properites  imply  that  the  vector  fields  Si,..., St.  denned  by 

(59)  s,=  Y,  sr 

are  mutually  orthonormal  vector  fields  in  V’F(i?"/27r^",(r“)  which  are  or¬ 
thogonal  to  Sq.  We  summarize  this  last  discussion  as  follows: 

Proposition  4  The  procedure  outlined  above  gives  rise  to  mutually  orthonc^- 
mal  vector  fields  Sii...,Sb  in  yF(il”/27r^",<r“)  which  are  orthogonal  to  Sq. 
The  corresponding  functions  defined  by  (44)  a  full  collection 

of  orthonormal  wavelets. 

3.6  Remarks 

Note  that  in  the  case  a  =  2  there  is  in  some  sense  a  “canonical”  method  of 
obtaining  Si  from  Sq.  This  is  essentially  the  formula  given  by  in  Subsection 
3.5.1.  Unfortunately,  as  is  quite  evident  from  the  construction  outlined  above, 
no  such  canonical  method  is  available  in  the  case  a  >  2. 

3.6.1  Cases  which  simplify 

In  certain  cases  the  general  procedure  outlined  above  is  not  only  unneces¬ 
sarily  cumbersome  but  may  not  produce  wavelets  with  desired  properties. 


45 


Sometimes  specialized  methods  may  be  simpler  and  produce  better  results. 
A  good  example  of  this  is  illustrated  by  the  method  of  obtaining  the  wavelets 
in  Subsection  3.7.1. 

Other  examples  arise  if  we  recall  that  the  second  step  of  the  recipe  is 
equivalent  to  finding  functions  . . . ,  Sb{0  which  are  periodic  and 

such  that 


(60) 


/  5o(e  -  27rB-^UQ)  SoiC  - 
Si(^  -  27r5- Vo)  5i(e  -  27rB-Vi) 


V  Sb{^  -  27rS-Vo)  Sb{^  -  27rB-Vi) 


5o(e-27rB-Vfc)  \ 
SiiC  -  27rB-^Ub) 

Sb{^-2TrB-'^Vb) 


is  an  orthogonal  matrix  for  all  ^  where  )Cb  —  {i^o,  •  •  •  >  *^6}  is  a  full  collection  of 
representatives  of  In  certain  cases  it  is  possible  to  do  this  directly. 

For  instance  if  a  =  2,  or  equivalently  6=1,  then  the  choice 

(61)  Si(«)  =  5o(«  - 


does  the  job,  where  mui  is  the  non-zero  element  in  K,a  and  k  is  any  any 
conveniently  chosen  element  of  This  is  equivalent  to  the  choice  of  Si 
given  in  Subsection  3.5.1.  Formula  (61)  is  essentially  the  classical  formula 
used  to  construct  univariate  wavelets  in  the  dyadic  case,  see  [14,  31,  33]. 

Another  such  example  arises  in  the  case  n  =  2,  A  =  B  =  21  when  50(0 
is  real  valued.  If 


then  one  can  verify  that 

Si(0  =  5o({  -  . 

52(0  =  S„((  -  , 

and 

53(0  =  5'o(e-7ri/3)e-‘>>’«> 

do  the  the  job,  see  [20,  33].  More  generally,  if  a  =  4  and  5o(0  is  real 
valued  an  analogous  choice  of  5i,  S^-,  a^nd  is  possible  whenever  2Z'^  J B2Z 
is  isomorphic  to  j 2  We  remind  the  reader  that  in  the  case  a  =  4  the 
group  Z'"  jBZ'^  can  be  isomorphic  to  one  of  two  abelian  groups,  .S'^/2.^^ 
or  ZjAZ. 


46 


3.6.2  Decay  properties 

A  desirable  feature  for  the  wavelets  is  that  they  enjoy  the  same 

decay  properties  as  the  scaling  function  which  they  are  associated  with.  For 
example  if  the  scaling  function  (j)  has  compact  support  or  exponential  decay 
then  it  would  be  nice  if  each  of  the  V’i’s  enjoyed  the  same  property.  Un¬ 
fortunately  the  general  recipe  outlined  above  does  not  necessarily  give  rise 
to  such  wavelets.  The  reason  for  this  is  that,  roughly  speaking,  if  a  >  2 
the  functions  ^i, . . . ,  ^6  produced  by  this  method  do  not  in  general  have  the 
same  smoothness  properties  as  Sq. 

There  are  many  specific  cases  where  one  can  directly  construct  wavelets 
which  enjoy  this  feature.  For  example,  it  is  not  difficult  to  find  compactly 
supported  wavelets  in  that  case  where  the  scaling  function  is  the  character¬ 
istic  function  of  a  compact  set,  see  Subsection  3.7.1.  On  the  other  hand  a 
comprehensive  treatment  of  this  issue  is  beyond  the  scope  of  this  presenta¬ 
tion.  Here  we  mention  several  fairly  wide  cases. 

The  case  a  =  2  This  case  is  well  understood  and  presents  no  difficulties. 
The  choice  for  given  by  (61)  gives  rise  to  a  wavelet  ip  =  tpi  with  es¬ 
sentially  the  same  decay  properties  as  the  corresponding  scaling  function  (p. 
For  example,  if  (p  has  compact  support  then  the  scaling  factor  Sq  must  be  a 
polynomial,  in  other  words  the  scaling  sequence  must  be  finite.  In  this  case 
the  function  given  by  (61)  is  also  a  polynomial  which  implies  that  the 
corresponding  wavelet  has  compact  support. 

The  case  a  =  3  The  general  recipe  outlined  in  Subsection  3.5.4  does  not  in 
general  give  rise  to  Si  and  S2  which  are  as  smooth  as  Sq.  On  the  other  hand 
the  procedure  indicated  in  Subsection  3.5.2  gives  rise  to  Si  and  S2  which 
have  the  following  properties: 

•  If  Sq  is  a  polynomial  or  is  analytic  then  and  S2  may  fail  to  have  this 
property.  However  both  and  S2  will  be  infinitely  differentiable. 

•  If  So  is  k  times  continously  differentiable  then  both  Si  and  S2  enjoy 
this  property. 
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The  case  a>  n  The  abstract  method  used  by  K.  Grochenig  [18],  see  also 
[33],  in  the  case  A  =  21  can  be  adapted  in  this  case  to  prove  the  existence  of 
Si,...,Sb  which  have  the  following  properties: 

•  If  So  is  a  polynomial  the  functions  Si,...,Sb  may  fail  to  be  polynomials. 
However,  they  will  be  analytic. 

•  If  50  is  analytic  or  k  times  continously  differentiable  then  functions 
5i, . . . ,  5b  enjoy  the  same  property. 

A  constructive  variant  of  this  method  has  been  derived  by  P.  Vial. 

Compact  support  The  desirability  of  nice  decay  properties,  particulary 
compact  support,  has  lead  in  the  univariate  dyadic  case  to  the  study  of 
“biorthogonal”  scaling  functions  and  wavelets,  see  [7, 42],  and  in  the  mutivari- 
ate  dyadic  case  to  a  detailed  study  of  various  pseudo  or  pre  scaling  functions 
and  wavelets  associated  with  box  splines,  see  for  example  [22,  25,  35,  38]. 


3.7  Examples 

3.7.1  Wavelets  generated  by  self  similar  sets 

Consider  the  multiresolution  analyses  discussed  in  Subsection  2.5.1.  Namely, 
suppose  ICa  =  {«o, . . . ,  Kb}  is  a  full  collection  of  representatives  of 
and  Q  is  the  compact  set  defined  by  (16)  which  is 

OO 

Q  =  {x  e  in'"  :x  =  Y^A~^€j,  €j  G  ICa)  • 

Furthermore  suppose  |Q|  =  1.  Then  xq  is  a  scaling  function  for  the  mul¬ 
tiresolution  analysis  it  generates.  In  this  case  a  fundamental  set  of  wavelets 
can  be  constructed  directly  as  follows: 

Recall  that  j  det  A\  =  a  and  b  =  a- 1.  Select  an  ax  a  orthogonal  matrix 
of  the  form 

/I  1  •••  1  \ 

1  «1,0  «1,1  •  •  •  ■“!,!» 

y/a  :  :  : 

\  Ubfl  Ub,l  ■  ■  ■  Ubjb  / 
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and  for  j  =  1, . . . ,  6  set 


b 

(62)  ^i(^)  ~  ^  "I  ^m)  • 

m=0 

The  functions  . . . ,  ^>6  are  a  full  collection  of  fundamental  wavelets.  The 
corresponding  orthonormal  wavelet  basis  for  L^{1R^) 

{a""^'^ipj{A'^x  -  A;)}jg{i,...,6},mez,  kez” 

is  a  natural  generalization  of  the  classical  Haar  basis  for  L'^(IR). 

Univariate  examples  If  A  =  2  and  =  {0, 1}  then  Q  is  the  interval 
[0, 1]  and  -tp  =  'ipi  is  the  function  in  Subsection  1.2.1.  The  collection 

{2m/2^(2m^  _  k)}k,meZ 

is  the  classical  Haar  basis  for  L‘^{JR). 

More  generally  if  A  =  N,  N  >  2,  and  JCa  =  {0;  •  •  •  >  “  1}  Q 

the  interval  [0, 1],  6  =  A/"  -  1  and  the  wavelets  ■01, . . . ,  '06  are  the  functions 
described  in  Subsection  1.4.1.  All  the  cases  when  Q  is  an  interval  give  rise  to 
similar  fundamental  wavelets  which  are  step  functions.  On  the  other  hand 
if  Q  is  not  an  interval  then  the  corresponding  wavelets  are  simple  functions 
which  are  somewhat  cumbersome  to  portray,  the  bivariate  analogues  are 
easier  to  describe  via  approprate  plots. 

Biv2iriate  examples  If  A  =  2/  and  the  members  of  =  {«o,  •  •  • » are 
described  by  the  columns  of 

/  0  1  0  1  A 

0  0  1  1  J 

then  Q  is  the  unit  square  [0, 1]^.  Of  course  there  are  many  possible  choices 
for  the  coefficients  in  represntation  (62).  The  popular  choice 


'  «1,0 

Ul,l 

Ul,2 

«1,3  ^ 

ri 

1 

-1 

-1  \ 

«2,0 

U2,l 

U2,2 

U2,3 

= 

1 1 

-1 

1 

-1 

V  «3,0 

U3,l 

^3,2 

“3,3  ) 

-1 

-1 

1  y 
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Figure  7:  The  Haar  like  wavelets  ^2  and 


Figure  8:  A  representation  of  the  wavelet  ^3  associated  to  the  tile  corre¬ 
sponding  to  A  =  2/  and  ICa  described  by  the  columns  of  ?? 

gives  rise  to  a  bivariate  analogue  of  the  Haar  system  with  the  fundamental 
wavelets  ipi,  ip2,  and  plotted  in  Figure  7. 

On  the  other  hand  if  the  members  of  )Ca  are  described  by  the  columns  of 

/  0  1  0  3  \ 

\  0  0  1  1 ; 

then  corresponding  tile  Q  together  with  the  coefficients  given  above  also 
gives  rise  to  three  wavelets  'tpi,  rp2,  and  V’3-  The  gray  level  plot  of  V’3  is 
displayed  in  Figure  8. 

If  the  matrix  A  and  the  set  of  digits  ICa  are  the  same  as  those  which 
give  rise  to  the  so-called  twin  dragon  set  Q  then  corresponding  fundamental 
wavelet  has  values  1  and  -1.  Its  gray  level  plot  is  displayed  in  Figure  9. 

We  finish  our  discussion  of  this  class  of  examples  by  considering  the  matrix 


'T 


Figure  9:  A  represntation  of  the  wavelet  associated  with  the  twin  dragon  tile. 


Figure  10:  A  representation  of  the  wavelets  in  the  last  tile  example. 

which  has  determinant  3.  If  the  members  of  1C  a  are  described  by  the  columns 

0  0-1 
0  1  1 

then  corresponding  tile  Q  together  with  the  coefficients  Uij  given  by 

(  t<l,0  »I,1  «1,2  )  =  V3  f  J  “1  S') 

V  “2,0  M2,1  «2,2  /  V  V2  ^  V2  J 

give  rise  to  fundamental  wavelets  tpi  and  ^2  whose  gray  level  plots  are  dis¬ 
played  in  Figure  10. 
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3.7.2  Wavelets  generated  by  univariate  polynomial  splines 


For  any  positive  integer  k  consider  the  dyadic  multiresolution  analysis  of 
L^{1R)  which  has  the  function  (j)  defined  by 

as  a  scaling  function.  See  Subsection  2.5.2.  The  corresponding  scaling  factor 
is  given  by 


q((:'\  —  I  k  +  47rm|  1  ^  _  1 


—im^ 


SfcC 

where  the  coefficients  {s^}  a^re  real  because  S  is  real  and  even.  Hence  the 
wavelet  ip  =  ipi  resulting  from  the  recipe,  see  (61),  can  be  expressed  as 

(63)  !>({)  =  e-‘^I^S{(m  -  Jr)^(f ) 


or 

</>(x)  =  E  -  m)  . 

tn^Z 

Note  that  for  any  positive  integer  k  the  sum 


^  le  +  27rml  2*= 

meZ 


can  expressed  in  terms  of  elementary  trigonometric  functions.  For  fc  =  1  this 
is  the  well  known  identity 


^  1^  +  27rm|  2  = 

m£Z 


4 

sin^^/2 


and  thus  for  fc  =  2, 3, . . ., 

1  4 

^|e  +  27rm|  (2A;-l)!dp-2sin2e/2  ’ 

Thus  the  Fourier  transforms,  (63),  of  these  fundamental  wavelets  can  be 
expressed  in  terms  of  elementary  functions.  Other  formulas  related  to  this 
example  can  be  found  in  [5,  14,  33]. 

Before  leaving  this  section  we  mention  that,  like  the  classical  Haar  system, 
these  wavelets  were  discovered  without  the  framework  provided  by  multires¬ 
olution  analysis.  They  were  discovered  by  J.  Stromberg  [40]  shortly  before 
this  framework  was  introduced. 


52 


3.7.3  Multivariate  analogues  of  the  univariate  spline  examples 

Here  we  indicate  some  of  the  details  only  for  the  polyharmonic  splines.  Var¬ 
ious  details  concerning  wavelets  arising  from  the  multiresolution  analyses 
generated  by  convolutional  iterates  of  characteristic  functions  can  be  found 
in  [39].  For  detailed  constructions  of  wavelets  and  psuedo  or  pre  wavelets 
generated  by  the  box  spline  multiresolution  analyses  see  [22,  25,  35,  38). 


Poly  harmonic  splines  For  any  integer  k  greater  that  n/4  consider  the 
multiresolution  analysis  of  associated  to  {2Z'^,A)  which  has  the  func¬ 

tion  6  defined  by 

=  -_JSL=_ 

\jY,iez''  IC  +  27rZ| 

as  a  scaling  function.  Assume  that  A  is  a  similarity  and  an  acceptable  dilation 
for  ZZ'^.  See  Subsection  2.5.3.  The  corresponding  scaling  factor  is 


(64) 


“  1  ie + 2™|-«  I 


where  B  =  A*. 

In  view  of  the  remark  recorded  immediately  after  the  statement  of  Propo¬ 
sition  3  the  fact  that  scaling  factor  (64)  is  real  valued  and  non-negative  im¬ 
plies  that  the  method  outlined  in  Subsection  3.5.3  can  be  used  to  to  obtain 
explicit  formulas  for  the  Fourier  transforms  of  a  collection  of  fundamental 
wavelets.  Like  the  scaling  function  ^  these  wavelets  enjoy  exponential  decay 
as  I  a:  I  goes  to  infinity. 

Note  that  in  the  case  a  =  2  the  formula  for  the  wavelet  ■0  =  ■0i  is 
analogous  to  formula  given  in  the  univariate  case. 


3.7.4  Compactly  supported  wavelets 

To  see  compactly  supported  wavelets  which  are  significantly  smoother  than 
the  examples  found  in  Subsection  3.7.1  consider  the  univariate  examples  cor¬ 
responding  to  the  scaling  functions  mentioned  in  Subsection  2.5.4.  The  sim¬ 
plest  specific  example  are  the  functions  0  and  0  described  in  Subsubsection 
1.2.3  and  plotted  in  Figure  2.  Since  a  =  2  in  these  cases,  an  application  of 
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(61)  gives  rise  to  a  wavelet  with  the  desired  properties;  that  is,  if 

2m+l 

Sj<f){2x  -  j) 

j=0 


then 


2m+l 

=  S  (-l)^S2m+i-i^(2a;  -  j) 
3=0 


where  the  k  in  formula  (61)  is  chosen  so  that  ip  has  support  in  the  same 
interval  as  <p. 

As  mentioned  in  Subsection  2.5.4  multivariate  examples  can  be  con¬ 
structed  by  taking  tensor  products  of  univariate  examples.  For  instance,  if 
(pi  and  ipi  are  scaling  function  and  wavelet  associated  to  a  univariate  dyadic 
multiresolution  anslysis  V^,  i  =  1,2,  then 

<^o(a^,y)  =  M^)<h{y) 


is  a  scaling  function  for  a  bivariate  dyadic  multiresolution  analysis  V  and 

‘Piix,y)  =  <Pi{^)^2iy) 

H>2{x,y)  =  'tpi{x)(p2{y) 

V’z{x,y)  =  ipi{x)ip2{y) 

are  a  full  collection  of  fundamental  wavelets  associated  with  V. 

As  indicated  in  Subsection  2.5.4  a  theory  of  smooth  compactly  supported 
wavelets  associated  with  a  general  dilation  matrix  A  has  not  yet  been  devel¬ 
oped.  We  close  this  section  with  the  following  amusing  example: 

Consider  the  acceptable  dilation  of 


and  the  bivariate  two  scale  difference  equation 

3 

¥^o(a^)  =  XI  Sj(po(Ax  -  kj) 
3=0 
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where  {soi  si,S2,  S3}  is  the  same  as  that  listed  in  Subsection  1.2.3  and 
{fco,  fci,  ^2,  fcs}  axe  the  columns  of 

/  0  1  2  2  \ 

\o  0  0  0  ;  ■ 

The  compactly  supported  solution  (po  of  this  difference  equation  is  a  scaling 
function  for  a  multiresolution  analysis  associated  to  {2Z^,A).  Both  (pQ  and 
an  associated  fundamental  wavelet 

3 

i=o 

can  be  expressed  in  terms  of  the  univariate  scaling  function  and  wavelet,  (j) 
and  V’)  decribed  in  Subsection  1.2.3,  namely 

(po{x)  = 

and 

V?i(rr)  =  ^(a;i)<^(x2) 

where  Xi  and  X2  are  the  first  and  second  conponents  of  x  respectively. 

3.7.5  Band  limited  wavelets 

Let  be  a  compact  subset  of  iR”  which  has  the  four  properties  listed  in 
the  beginning  of  Subsection  2.5.5  and  let  V  be  the  multiresolution  analysis 
associated  to  A)  which  has  where  0(0  =  as  a  scaling  function. 

Suppose  are  mutually  essentially  disjoint  subsets  of  BQ  each  of 

which  which  are  congruent  to  0  mod  2Tr2Z'^  and  essentially  disjoint  from  fl. 
In  other  words,  if  Oq  =  ^  then 

•  flj  C  BQ,  for  each  i  G  {0, 1, . . . ,  6}. 

•  +  27rA:}|  =  0  for  any  element  k  in  \  {0}  and  each  i  G 
{0,1,..., 6}. 

•  For  each  i  G  (0, 1, . . . ,  6} 

U  {fii  +  2Trk}  =  nr . 

jfceZ" 
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Finally,  for  each  i  e  {1, . . . ,  consider  the  function  ^  defined  by 

Mo  =  xcliiO  ■ 

Then  the  functions  , . . . ,  are  a  collection  of  fundamental  wavelets  asso¬ 
ciated  with  V. 

Just  as  the  scaling  function  ^  these  wavelets  have  poor  decay  properties 
at  infinity. 

The  univariate  case  where  A  =  B  =  2  and  <f){x)  =  sincx  the  correspond¬ 
ing  wavelet  tpi  is  the  sine  wavelet  mentioned  in  Subsection  1.2.2  and  plotted 
in  Figure  1. 

In  the  bivariate  case  where 


and  4>{x)  =  sincx  then  the  corresponding  wavelet  ipi  is  the  sine  wavelet  ip 
mentioned  in  Subsection  1.4.2  and  plotted  in  Figure  3. 

Wavelets  with  better  decay  properties  These,  of  course,  are  derived 
from  the  scaling  functions  with  better  decay  properties  mentioned  in  Subsec¬ 
tion  2.5.5.  Unfortunately  their  construction  is  not  as  simple  as  that  found 
above.  In  the  case  a  =  6-|-l  =  2we  can  apply  the  canonical  procedure  which 
results  from  the  use  of  formula  (61).  Fortunely,  since  the  scaling  factor  is 
non-negative,  in  the  case  o  >  2  the  procedure  outlined  in  Subsection  3.5.3 
can  be  applied  to  obtain  the  SPs  to  get  the  corresponding  ipPs.  In  either  case 
the  wavelets  will  enjoy  the  same  decay  properties  as  the  scaling  function. 
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Abstract 

It  is  well  known  that  the  integer  translates  of  the  scaling  func¬ 
tion  associated  with  a  given  scaling  sequence  may  fail  to  be  mutually 
orthogonal.  Here,  we  address  seversd  technical  questions  related  to 
this  phenomenon.  For  instance,  we  show  that  the  scaling  function 
naturally  associated  with  a  finite  scaling  sequence  always  generates  a 
mxdtiresolution  analysis  and  give  elementary  but  non-trivial  examples 
of  scaling  sequences  which  give  rise  to  pathological  scaling  functions. 


1  Introduction 


Suppose  if  satisfies  the  two  scale  difference  equation 

(1)  ¥’(®)  =  X)  -  fc), 

k 


where  x  is  any  real  number  and  {sjt}  is  a  finite  sequence  of  scalars.  Then  a 
necessary  condition  that  {<p{x  —  ky}kez  ^  orthonormal  system  of  func¬ 
tions,  namely 

(2) 


/OO  - 

(p{x  —  m)(p{x  —  n)dx  =  Sm,r 

-OO 


which  generates  a  dyadic  multiresolution  analysis  is  that  the  sequence  {sfc} 
satisfies  the  following  constraints; 


(3) 

X)  ^kSk-2n  —  2^0, n 
k 

and 

(4) 

XI  Sfc  =  2, 

1 


where  8m, n  is  the  Kronecker  delta.  The  sequence  {^jt}  is  often  referred  to  as 
a  scaling  sequence  and,  if  both  (1)  and  (2)  hold,  the  function  <p  is  called  a 
scaling  function.  Relation  (3)  is  a  consequence  of  (1)  and  (2)  while  relation 
(4)  follows  from 

rcxj 

(p{x)dx  =  1, 


/: 


a  condition  which  is  imphed  whenever  the  closed  linear  span  of  {^{x  —  k)}k^z 
generates  a  dyadic  multiresolution  analysis.  See  [4]  for  more  details. 


1.1 

Recall  that  if  {s*,}  is  a  finite  scaling  sequence  which  satisfies  (3)  and  (4)  the 
canonical  solution  of  (1)  is  given  by  the  function  tp  which  is  well  defined  by 
the  formula  for  its  Fourier  transform 


(5) 

where 

(6) 

and 

(7) 


■fd)  =  n 

k=l 

^  k 

<P{0  =  / 

•/— OO 


Unfortunately,  examples  show  that  (p  need  not  satisfy  (2).  See  Subsection 


1.4. 

On  the  other  hand  such  a  p  always  generates  a  dyadic  multiresolution 
analysis.  To  be  more  precise  let  V{}p)  be  the  closed  hnear  span  of  {y>(s  — 
k)}k&z  in  L^{IR)  ■  We  say  that  p  generates  a  multiresolution  analysis  when¬ 
ever  there  is  a  multiresolution  analysis 

V  =  {...  C  U_i  C  Uo  C  C  Uz  C  ...} 

with  Vo  =  y{v)- 

Theorem  1  Suppose  {s^}  is  a  finite  scaling  sequence  which  satisfies  ( 3)  and 
(4)  andp  is  the  function  which  is  described  by  (5).  Then  p  generates  a  dyadic 
multiresolution  analysis. 

If  p  fails  to  satisfy  (2)  then  it  is  somewhat  technical  and  awkward  to 
describe  a  scaling  function  and  corresponding  scaling  coefficients  for  the  mul¬ 
tiresolution  analysis  which  is  generated  by  p.  The  Fourier  transform  of  such 
a  scaling  function  po  is  essentially  given  by 


<fo{0  = 


TiO 
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I£  all  the  coefficients  {s*,}  axe  real,  this  multiresolution  analysis  always  has 
a  scaling  function  which  is  supported  in  the  interval  [0,  oo)  ;  it  also  has  one 
supported  in  (— oo,0].  The  proof  of  the  Theorem  and  the  details  to  these 
remaxks  can  be  found  in  Subsection  2.1. 

1.2 

The  conclusion  of  Theorem  1  may  fail  if  {sfc}  is  not  assumed  to  be  fimte.  For 
example,  there  axe  sequences  {s*,}  which  satisfy  both  (3)  and  (4)  and  decay 
faster  than  the  reciprocal  of  any  polynomial  such  that  the  corresponding 
function  y  given  by  (5)  fails  to  generate  a  multiresolution  analysis. 

Nevertheless  this  theorem  can  be  extended  to  sequences  which  axe  not 
finite  and/or  do  not  necessarily  satisfy  (3).  For  example  the  arguments  in 
Subsection  2.1  show  that  the  conclusion  of  Theorem  1  holds  whenever  the  bi- 
infinite  sequence  {sk}  satisfies  (4)  and  enjoys  exponential  decay  as  \k\  — >  oo. 

1.3 

If  {^fc}  is  a  scaling  sequence,  consider  the  corresponding  periodic  function 
S{C)  defined  by  (6).  Recall  that  the  constraints  (3)  and  (4)  on  {sk}  are 
eqrdvalent  to  the  conditions 

(8)  i5(oi^+i5'(^+^)r  =  i 

and 

(9)  S(0)  =  1 
on  S{^).  Furthermore,  relation  (1)  becomes 

(10)  m  =  simncm, 

and  relation  (2)  is  equivalent  to 

OO 

(11)  = 

fc=— oo 

The  interplay  between  {sfc},  5(^),  and  the  scaling  equation  (1)  may  be  made 
more  succinct  if  to  {«*}  one  associates  the  train  of  Dirac  delta  distributions 


(12) 

k 

Then 

(13) 

^ 

3 


and  equations  (1)  and  (10)  become 

(14)  (p{x)  =  s*  (p{2x) 

and 

(15) 


respectively. 

Given  a  scaling  sequence  {sjs;}  which  satisfies  the  constraints  (3)  and  (4) 
and  the  corresponding  periodic  function  5'(^)  defined  by  (6)  consider  the 
function  Suiji)  defined  by 


(16) 


Sn((.)  =  S{Nt), 


where  is  a  positive  integer.  Clearly,  5//  is  also  a  2ir  periodic  function  and 
if  N  is  odd  it  also  satisfies  (8)  and  (9).  In  other  words 


(17) 


iMor+iMe+^)p  =  i 

5;(r(0)  =  1 


and 
(18) 

whenever  iST  is  an  odd  integer. 

If  s(s)  is  the  train  of  delta  functions  which  corresponds  to  {s*}  and  5(^) 
then  the  train  sn{x)  corresponding  to  Sn{0  given  by 


(19) 


sn{x)  =  J^six/N), 


which,  in  view  of  the  fact  that  S{x/N)  =  NS{x),  may  be  expressed  as 
(20)  sn{x)  =  '^SkS{x  —  Nk). 

k 

The  corresponding  scaling  sequence  may  be  somewhat  awkwardly 

described  by 


(21) 

In  other  words 
(22) 


_  j  Sfn  if  fe  =  =  0,  ±1,  ±2, ... 

1  0  otherwise. 


In  view  of  the  relationship  (21)  between  {5*}  and  {5j\r,fc}fc  1  we  say  that 
{siv^,fc}fc  is  a  dilate  of  {sfc}.  Similarly,  we  say  that  5jv  is  a  dilate  of  S. 

If  (p{x)  is  the  canonical  scaling  function  given  by  (5)  associated  with  {sk} 
then  the  function  defined  by 


(23) 


4 


is  the  corresponding  scaling  function  associated  with  the  sequence 

If  JV  is  an  odd  integer,  in  view  of  (17)  and  (18),  the  sequence  {sN,k}k  is 
also  a  scaling  sequence  which  satisfies  (3)  and  (4).  Unfortunately  if  JV  >  1, 
the  corresponding  scaling  function  (pjf  necessarily  fails  to  satisfy  (2),  namely, 
the  collection  {^jv’(®  —  ^)}fce^  fails  to  be  orthogonal.  See  Section  2.3  for  more 
details. 

The  above  observations  give  rise  to  an  easy  method  for  constructing  scal¬ 
ing  sequences  which  satisfy  (3)  and  (4)  but  whose  corresponding  scaling  func¬ 
tions  given  by  (5)  fail  to  satisfy  (2). 

1.4 

In  the  previous  section  we  described  a  method  for  constructing  scaling  se¬ 
quences  which 

•  satisfy  (3)  and  (4) 

■  fail  to  give  rise  via  (5)  to  scaling  functions  (p  which  satisfies  (2). 

It  is  not  immediately  clear  whether  there  are  scaling  sequences  which 
enjoy  these  properties  and  which  do  not  arise  in  this  way.  In  view  of  this  it 
is  natural  to  ask  whether  every  sequence  which  enjoys  the  properties  listed 
above  is  an  appropriate  dilate  {5//,fe},  with  N  odd  and  >  1,  of  another  scaling 
sequence  {sfc}. 

This  question  was  answered  in  the  negative  by  Cohen  and  Sun  [3].  Be¬ 
cause  their  construction  also  resolves  a  more  delicate  question,  see  Subsection 
1.5,  it  is  unnecessarily  compUcated.  Here  we  present  an  alternate,  more  ele¬ 
mentary  solution. 

Observe  that  the  existence  of  a  finite  scaling  sequence  {s^}  which  has  the 
appropriate  properties  and  is  not  a  dilate  of  another  scaling  sequence  follows 
from  the  existence  of  a  trigonometric  polynomial 

(24)  S(«  =  E 

fc=0 

which  satisfies  (8)  and  (9),  gives  rise  via  (5)  to  a  scaling  function  p  which 
fails  to  satisfy  (2),  and  has  the  property  that  S{C/N)  is  not  2tv  periodic  for 
any  integer  N  >  1. 

Proposition  1  For  every  even  integer  n  which  is  >  8  there  are  trigonomet¬ 
ric  polynomials  S  of  the  form  (24)  which  enjoy  the  following  properties: 

(i)  The  coefficients  {sfc}  are  real,  sq  ^  0,  and  s„_i  ^  0. 

(ii)  S{^)  satisfies  (8)  and  (9). 
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(iii)  The  function  (p  defined  via  (5)  fails  to  satisfy  (2). 

(iv)  S{^/N)  is  not  2t:  periodic  for  any  integer  N  >  1. 

Our  construction  of  the  polynomials  S  whose  existence  is  guaranteed  by 
the  above  proposition  is  based  on  the  following  lemma 

Lemma  1  If  S  is  a  2ir  periodic  trigonometric  polynomial  and  has  1  + 
as  a  factor,  where  N  is  an  odd  integer  greater  than  1,  then  the  function  (p 
defined  hy  (5)  fails  to  satisfy  (2). 

The  idea  is  to  construct  a  2tv  periodic  trigonometric  polynomial  S  which 
satisfies  (8),  (9),  and  the  conditions  of  Lemma  1  but  is  not  a  dilate  of  another 
27r  periodic  trigonometric  polynomial.  In  the  case  N  =  Z,  take 

(28)  S(()  = 

where  F{^)  is  a  trigonometric  polynomial  chosen  so  that  S  satisfies  (8)  and 
(9),  Of  course  the  choice  F{()  —  1  does  the  job  but  the  resulting  5(^)  is  a 
dilate  of  — . 

Lemma  2  For  every  even  integer  n  which  is  >  4  there  are  trigonometric 
polynomials  F(^)  of  the  form 

(26)  F(i)  =  E 

k—0 


such  that 

•  cq  ^  0  and  Cn  ^  0 

m  The  polynomial  S  defined  by  (25)  satisfies  (8)  and  (9) 

•  This  polynomial  can  he  chosen  so  that  Ci  ^  0  . 

The  last  item  of  the  Lemma  2  guarantees  that  S{^fN)  is  not  27r  periodic 
for  any  integer  N  >  1. 

The  coefficients  {c*,}  in  (26)  axe  solutions  to  a  system  of  quadratic  equa¬ 
tions.  A  specific  numerical  solution  to  this  system  in  the  case  n  =  4  is  given 

ty 

{co,  Cl,  C2,  C3,  C4}  =  {0.280046,  -0.553373, 0.609149, 0.441001, 0.223178}. 

This  gives  rise  to  a  polynomial  <S'(^)  in  the  form  (24)  with  n  =  8  and  which 
satisfies  the  conclusion  of  Proposition  1.  The  coefficients  {sq,  ...,57}  3-re 

{0.140023,  -0.276687, 0.304574, 0.360523, 

-0.165098, 0.304574, 0.220500, 0.111589}. 


Details  are  in  Section  2.4  . 
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1.5 


In  view  of  Lemma  1  one  may  a^k  whether  the  presence  of  the  factor  1  + 
N  odd  and  greater  than  1,  is  necessary  in  the  polynomials  S{^)  whose 
existence  is  guaranteed  by  Proposition  1.  The  example  constructed  by  Cohen 
and  Sun  [3]  shows  that  the  answer  is  no. 

Here,  by  slightly  modifying  the  example  given  in  [3],  we  construct  the 
polynomial  with  real  coefficients  of  minimal  degree  which  does  not  have  the 
factor  1  +  N  odd  and  greater  than  1,  and  enjoys  properties  (i)-(iv)  of 

Proposition  1. 

Consider  the  2x  periodic  trigonometric  polynomial 

+  1 
2 


whose  roots  in  (— x,7r]  are 


— Ttt  — Stt  — 37r  — TT  TT  37r  Stt  Ttt  . 
^  9  ’  'Y'  9’  T’  T’ 


This  polynomial  has  factors 


and 


G2{0  = 


GjQ 

Gi(0  ■ 


Observe  that  the  set  of  roots  of  <j2(^)  iii  is  the  set 


fi  =  { 


— Ttt 


—Sir  — TT  TT  Stt  77r 


Define  the  polynomial  H{()  via 


(27)  H{0  =  |G2(0r  {1^1(01^  + c|G2(^  +  7r)P  cose} 


where  the  real  constant  c  is  chosen  so  that  the  expression  in  braces  is  positive. 
This  polynomial  enjoys  the  following  properties: 

•  i?(e)  is  a  non-negative  27r  periodic  trigonometric  polynomial  whose  set 
of  roots  in  (— tt,  tt]  is  the  set  D. 

.  H{-()  =  H{0 . 

.  H(0  +  /fW  +  >r)=l. 
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The  last  two  items  are  can  be  easily  seen  by  first  re-expressing  H{()  as 

m)  =  iG(or + c|G2(opig2(^ + cos^ 

and  observing  the  behavior  of  each  of  the  two  summands.  These  properties 
imply  that 

15 

H{^)  =  ^ak  cos  kC 
k-O 

which,  in  view  of  the  lemma  of  Riesz,  can  be  factored  as 


(28) 

where 

(29) 

with  real  coefficients  Sk- 


m)  =  i-s-cop 


15 


S«)  =  E  «<= 

fc=0 


-iH 


Proposition  2  The  27r  periodic  polynomial  S{^)  defined  by  (28)  and  (29) 
satisfies  properties  (i)-(iv)  of  Proposition  1  and  does  not  have  1  -|-  e  as  a 
factor,  for  any  integer  N  >  1. 


1.6 

Up  to  this  point  we  have  been  mainly  concerned  with  the  canonical  solution 
of  (1),  that  is,  the  solution  (p  described  by  (5).  Since  there  are  other  solutions 
of  (1),  which  are  not  constant  multiples  of  the  canonical  solution,  it  is  natural 
to  ask  when  such  solutions  satisfy  (2).  The  following  theorem  should  resolve 
most  of  these  questions. 

Theorem  2  Suppose  {s*}  is  a  finite  scaling  sequence  which  satisfies  (3)  and 
(4)  and  (f  is  the  function  whose  formula  is  given  by  (5). 

•  If  ip  fails  to  satisfy  (2)  then  every  solution  of  (1)  also  fails  to  satisfy 

(2). 

•  If  ip  satisfies  (2)  then  X  is  another  solution  of  (1)  which  satisfies  (2)  if 
and  only  if 

m  =  mm 

where  h{2^)  =  h{^)  and  |/i(^)|  =  1  almost  everywhere  on  M. 

This  result  may  be  roughly  rephrased  as  follows:  If  X  is  a  solution  of  (1) 
which  satisfies  (2)  then  the  canonical  solution  p  also  satisfies  (2)  and 

|A(01  =  1^(01 

almost  everywhere  on  St. 
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1.7 


Finally  we  mention  that  the  examples  in  Subsections  1.4  and  1.5  have  lowest 
degree  possible.  More  precisely, 

Proposition  3  Every  trigonometric  polynomial  S  of  the  form  (24)  which 
satisfies  the  conditions  of  Proposition  1  has  n  >  8.  If,  in  addition,  S  is 
required  to  not  have  1  +  N  odd  and  greater  than  1,  as  a  factor  then 

n  >  16. 


2  Details 

These  sections  contain  details  omitted  in  the  previous  sections.  We  use  the 
convention,  employed  in  [1],  that  the  details  to  Section  l.n  are  contained  in 
Section  2.n  . 


2.1 

We  begin  by  recalling  the  following  facts  from  [7]  : 

Proposition  4  A  function  (p  in  L\M)  generates  a  dyadic  multiresolution 
analysis  if  and  only  if  (p  enjoys  the  following  properties  : 

•  There  is  a  2'X  periodic  measurable  function  S{^)  such  that 

(SO)  m  =  s((/2)vm)- 

•  For  almost  all  ^ 

OO 

(31)  E  |^(C-27rA!)P>0. 

fe=:— OO 

•  For  every  finite  interval  I 

(Sg)  Jim  /  { 1^(01“  /  E  -  2x^)1%  =  1. 

J-oo  j2  3I  \  J2-SI  t  f 

where  2~^I  =  {x  :  2^x  €  /}  . 

Proposition  5  Suppose  <p  is  a  function  in  L^{JR)  which  enjoys  the  properties 
listed  in  Proposition  4  ("■nd  P(^)  is  a  measurable  27r  periodic  function  such 
that 

OO 

i/’(e)p=  E  iv«-2x^)r 

k=—oo 

Then  the  function  po  defined  by 

(SS)  MO  = 

is  a  scaling  function  for  the  multiresolution  analysis  generated  by  p. 
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Suppose  {sk}  and  <p  are  as  in  the  hypothesis  of  Theorem  1.  Then  9?  is  in 
L^[1R)^  has  compact  support  and  satisfies  (30)  with  S  given  by  (6);  see  [4]. 
Define  the  periodic  function  $  via 


00 

(34)  m=  E  l^(£-2rf)|’. 

k——oo 


By  virtue  of  the  Poisson  summation  formula,  $  is  a  trigonometric  polynomial 
and  hence  $(^)  >  0  almost  everywhere. 

Since  S  satisfies  (8)  and  (9),  it  follows  that  S{Tr)  =  0  and  hence,  in  view 
of  (5),  yj(27rfc)  =  0  for  all  A:  €  ^\{0}.  Thus  1^(0IV^(0  is  continuous  in  a 
neighborhood  of  zero  and  equals  one  for  ^  =  0.  This  imphes  that  (32)  holds 
and,  together  with  the  previous  observations,  allows  us  to  conclude  that  all 
the  conditions  in  Proposition  4  are  satisfied,  so  the  proof  of  Theorem  1  is 
complete. 

Observe  that 

.  m  saimm  simpimMtm 

MO  - 

so  that  (fa  satisfies  the  scaling  equation 

(35)  MO  =  ■yo(£/2)^o(£/2), 


where 

(36) 


SoiO  = 


SiQPiO 
P{20  • 


Note  that  So{C)  satisfies  (8)  and  (9)  as  it  should.  Also  note  that  if  ip  satisfies 
(2)  then  $(^)  =  1  so  that  P  can  also  be  taken  to  be  =  1  and,  in  this  case, 

V^o  =  P-  - 

For  the  function  P  in  (33)  one  can  simply  take  P{i)  =  ^ $(^),  the  positive 
square  root  of  $(|).  If  $(0  =  ^(— example  this  is  the  case  when  all 
the  coefficients  {sk}  are  real,  then  by  virtue  of  the  Riesz  lemma,  see  [4],  P(^) 
can  be  chosen  to  be  a  trigonometric  polynomial;  in  this  case  cr(e”*^)  =  5'o(^) 
is  a  rational  function  of  z  =  e“*^. 

Let  P(^)  =  p(e“*^)  be  a  27r  periodic  trigonometric  polynomial  satisfying 


(37)  |P(0P  =  ^(0 

and  observe  that  the  following  are  also  27r  periodic  trigonometric  polynomials 
satisfying  (37). 

•  e*”^P(^)  for  any  integer  m  and 

•  P(^)(l  —  ae“*^)/(e“*^  —  a)  where  o  is  a  root  of  p(z). 
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The  last  item  impUes  that  P{^)  =  p(6  can  be  chosen  so  that  aU  the  roots 
a  oip{z)  satisfy  either  ja]  <  1  or  jaj  >  1. 

If  P  is  chosen  so  that  all  the  roots  of  p{z)  are  >  1  then,  in  view  of  the  fact 
that  \So{0\  <  aU  the  poles  b  of  the  rational  function  cr(z),  if  any,  satisfy 
|61  >  1  and,  by  multiplying  P(^)  by  additional  factors  of  if  necessary,  we 
see  that  5'o(^)  can  be  expressed  as 

OO 

5o(«  = 

ik=0 

where  hmsup;i,_^oo  <  1*  In  tins  case,  since  tpo{x)  is  the  solution  of  the 

two  scale  difference  equation 

OO 

<Po{x)  =  ^*ipo(2a;  -  k) 

k-O 

via  fixed-point  iteration  starting  with  the  indicator  function  of  the  interval 
[0, 1],  it  follows  that  <po{x)  has  support  in  [0,  oo). 

An  analogous  computation  produces  y>o(®)  with  support  in  (— oo,0]  . 


2.2 

To  see  that  the  conclusion  of  Theorem  1  may  fail  if  {sfc}  does  not  decay 
sufficiently  rapidly,  consider  the  27r  periodic  function  S{^)  which  is  defined  as 
follows;  begin  with  any  non-negative  C°°  function,  h{^),  which  is  supported 
in  /e  =  ;  1^1  <  -H  c}  and  is  positive  on  1^/2',  here  e  satisfies  0  <  e  <  tt. 

Define  g{()  by 

^  MO 

and  note  that 

(38)  E  +  = 

k&Z 

Let 

^’(0  =  E 

kez 

and  finally  define  S{()  via 


S{0  =  T{2N0 


where  JV  is  a  positive  odd  integer.  Since  S{^)  is  an  infinitely  differentiable 
27r  periodic  function,  we  may  write 

(39)  S(f)  =  i  E 

^  kez 
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where  the  sequence  {sfc}  decays  faster  than  the  reciprocal  of  any  polynomial. 
In  other  words,  {sfc}  enjoys  the  estimates 

N<ap(i+ifcir 

for  all  positve  p  where  the  constant  Gp  is  independent  of  k.  In  view  of  (38) 
and  the  fact  that 

|rK)|=  =  E  l9(f  +  . 

kez 

the  function  S{()  satisfies  (8)  and  (9)  which  implies  that  {s^}  satisfies  (3) 
and  (4). 

Now,  if  <p  is  the  function  defined  in  terms  of  S{^)  via  (5)  then  ip  vanishes 
in  a  neighborhood  of  2fc7r/iV,  for  all  fc  in  .^  \  {0}.  Thus  if  fV  is  greater  than 
1  then,  recalling  that  N  is  odd,  it  follows  that 

l^(^  +  27rA:)l2 

kez 

vanishes  in  a  neighborhood  of  2x/JV.  In  view  of  the  second  item  in  Proposi¬ 
tion  4,  we  may  conclude  that  (p  fails  to  generate  a  multiresolution  analysis. 


2.3 

To  see  that  ipN(x)  defined  by  (23),  with  an  odd  integer  N  >  I,  necessarily 
fails  to  satisfy  (2),  assume  that  {sfc}  is  finite  and  observe  that  y>(27rm)  =  0 
for  all  m  in  \  {0}.  Thus  ^jv^(^)  =  <p{N^)  vanishes  whenever  C  =  2m-KlN 
for  zilm  'va.  Z\  {0}.  Hence 

OO 

k——oo 

is  a  continuous  periodic  function  which  vanishes  when  ^  =  {2Tr  fN)  -1-  27rm, 
m  E  2Z  and  the  desired  result  follows  since  the  orthogonality  of  the  system 
{^jv^(a:  —  k)}k  is  equivalent  to  $j\r(0  =constant. 

The  general  case  follows  from  a  similar  argument. 


2.4 

To  see  Lemma  1  write 

(40)  S(()  =  - F(0 

where  i^(0)  =  1  and  observe  that 


(41) 


<p{0  - 


1  - 

2iN$ 


HO 
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where 

OO 

m = n 

k=l 

is  an  analytic  function.  Prom  (41)  we  may  conclude  that  vanishes  when¬ 
ever  i  =  2rmrlN  for  all  m  in  \  {0}  and  the  desired  result  follows  from  the 
fact  that 

OO 

*«)=  E  iip«-2ri)r 

k——oo 

vanishes  when  ^  =  (2irfN)  27rm,  m  E  Z. 

To  find  a  suitable  polynomial  F  for  Lemma  2  substitute  (25)  into  (8)  to 
get 

(42)  iG(on^(oi^ + m + ^)i  =  1, 

where  16?(^)|^  =  (cos(^))^  and,  of  course,  |G(^  +  7r)p  =  (sin(^))^  Recall 
that 

|G(«P  +  |G(e+’f)P  =  i 

and  observe  that  this  together  with  (42)  imphes 

(43)  iG(^)r{i^’U)P  - 1} + m + ^)i^{ine + ^)p  - 1}  =  0.  ■ 

Now  (43)  is  not  difficult  to  solve  for  1^(0 P  if  we  choose  it  to  be  of  the  form 

(44)  \FiO?  =  \GU  +  ^WQiC  +  ^)  +  ^. 
where  Q  is  a  trigonometric  polynomial.  Thus 

m+^)\^  =  mwQio  +  T^ 

and  substituting  the  last  two  expressions  in  (43)  results  in 

1G(0PIG'(^  +  ^)P{QU  +  ^)  +  QiO}  =  0 

or,  more  simply, 

(45)  <3(^  +  7r)-l-Q(e)  =  0. 

It  is  cleax  that  any  linear  combination  of  cos((2fc  —  1)^),  fc  =  0, 1, 2, . . .,  solves 

(45) .  Namely  every  trigonometric  polynomial  Q  of  the  form 

i 

(46)  QiO  =  XI  cos((2A!  -  1)^) 

fc=i 

is  a  solution  of  (45).  In  particular  choosing 

(47)  g(^)  =  4acos(^) 
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results  in 

(48)  I'P’COr  —  o(cos4^  +  cos  2^  —  2cos^)  +  1. 

Details  of  this  calculation  can  be  found  at  the  end  of  this  section. 

In  general  it  is  clear  that  if  the  absolute  values  of  the  coefficients  in 
representation  (46)  are  sufficiently  small  then  the  right  hand  side  of  (44)  is 
positive.  In  this  case  the  Riesz  lemma  mentioned  above  implies  the  existence 
of  a  trigonometric  polynomial 

2f+2 

(49)  F{C)  =  ^ 

fc=0 

which  satisfies  (44).  In  particular  almost  all  choices  of  sufficiently  small 
in  representation  (46)  will  reult  in  non-zero  cq,  ci,  and  C21+2  in  (49).  In  such 
cases  the  polynomial  S{i)  described  by  (25)  satisfies  both  (3)  and  (4)  and 
S{^/N)  fails  to  be  27r  periodic  for  every  integer  N  ^  1. 

Returning  to  (48)  note  that  the  maximum  of  cos  4^  -|-  cos  2^  —  2cos^ 
occurs  at  ^  =  TT  and  is  4.  Numerical  methods  show  that  the  minimum  occurs 
at  ^  =  0.818919  and  is  -2.42404;  these  estimates  are  accurate  to  the  listed 
number  of  decimal  places.  Hence  the  right  hand  side  of  (48)  is  positive  if 

(50)  -  0.25  <  a  <  0.412534. 

For  such  a,  the  Riesz  lemma  [4,  page  172]  implies  the  existence  of  the  desired 
trigonometric  polynomial 

(51)  = 

which  satisfies  (48).  If  a  ^  0,  then  it  is  clear  that  both  coeficients  cq  and  C4 
in  (51)  are  not  equal  to  0  and  the  polynomial  S{C)  described  by  (24)  has  the 
property  that  S{^/N)  fails  to  be  27r  periodic  for  every  integer  N  ^  1. 

To  get  an  exphcit  collection  of  the  coefficients  c^  for  the  polynomial  F(^) 
given  by  (51)  observe  that  such  coefficients  must  satisfy  the  system  of  equa¬ 
tions 

Cq  C^  -i  Cj  -t-  Cg  4-  C4  =  1 

CqCi  -i  C1C2  "i  C2C3  -f  C3C4  =  —a 
C0C2  +  C1C3  +  C2C4  =  a/2 

C0C3  +  C1C4  =  0 

C0C4  =  a/2 

where  a  ^  0  and  satisfies  the  constraints  (50).  This  system  of  equations  can 
be  easily  solved  by  numerical  method  such  as  Newton’s  method.  In  the  case 
a  =  0.125  one  such  solution  is 

{co,  Cl,  C2,  C3,  C4}  =  {0.280046,  -0.553373, 0.609149, 0.441001, 0.223178}. 
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To  see  (48)  recall  that 


so  that  if  Q  is  given  by  (47) 

|G(0P<3(0  =  2a(cos^  +  cos  3^cos^). 


Now  using 

we  get 
cind 


cos  3^  cos  ^  =  —(cos  4^  +  cos  2^) 

1G(0I^^(0  =  a(cos4^  +  cos  2^  +  2cos^) 

+  ^)  =  a(cos4^  +  cos  2^  -  2cos^). 


Substituting  the  last  formula  into  (44)  gives  the  desired  result. 


2.5 

To  see  that  the  constant  c  can  be  chosen  so  that  the  expression  in  braces  in 

formula  (27)  is  positive,  note  that  |(ji(C)P  is  positive  on  (-7r,ir]  except  at 

4  =  7r/3  and  at  ^  =  — ?r/3.  Since  the  set  of  roots  of  |C?2(^  +  ^)1^  (—' is 

the  set  .  ^  ^  .  c 

^  ^  Stt  — 47r  — 27r  ^  27r  4x  8x., 

0,  y ,  y ,  y> 

there  is  a  positive  e  such  that  |G2(^  +  x)p  is  positive  for  ^  in  Jg  =  (x/3  — 
e,  x/3  +  e)  U  (5x/3  -  e,  5x/3  +  e)-  Now,  if 


m  =  inf 


\Giiono 


and 

M—  sup  |G2(^  +  x)p. 

46[0,2^]\Jg 

then  any  c  such  that  0  <  c  <  mfM  will  do  the  job. 

To  see  Proposition  2  note  that  it  follows  from  the  properties  of  H{^) 
that  S{^)  satisfies  the  first,  second  and  fourth  conditions  of  Proposition  1  . 
Furthermore,  since  fi  is  the  set  of  roots  of  5(^)  in  (— x,  x],  <S'(^)  does  not  have 
1  -1-  as  a  factor.  However,  the  third  condition  of  Proposition  1  needs 

some  explanation. 

Observe  that  the  set  fl  +  x  fisted  above  has  the  following  properties: 

•  +  x  is  a  non-trivial  invariant  cycle  in  (— x,  x]  for  the  operation  ^  i->  2^ 

mod  2x. 
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•  |i?(^)|  =  1  for  all  ^  in  ft  +  TT. 

In  view  of  these  properties  Cohen’s  result  concerning  such  cycles  implies  that 
S{i)  satisfies  the  third  condition  of  Proposition  1.  See  [4,  page  188]  or  [2,  3,  5] 
for  more  details. 

Before  closing  this  subsection  we  mention  that  the  construction  of  5(^) 
here  is  really  analogous  to  that  in  Subsection  2.4  with  |(?i(^)P  playing  the 
role  here  that  1  played  there. 


2.6 

To  see  Theorem  2,  suppose  is  a  finite  scaling  sequence  which  satisfies 
(3)  and  (4),  (p  is  the  function  whose  formula  is  given  by  (5),  and  A  is  another 
solution  to  (1).  Then 

(52)  m  =  mm 

where  h  is  a  function  which  satisfies 

(53)  K20  =  h{0 

almost  everywhere  on  M.  Observe,  that  if  A  is  essentially  bounded  then  so 
is  h,  namely, 

(54)  ess  sup  |fe(^)|  <  ess  sup|A(^)|. 

Inequality  (54)  follows  from 

ess  sup  1A(^)|  >  ess  sup  \h{$)\  mf  |^(^)| , 

|^|<5 


the  fact  that  infi^|<5  |^(0I  arbitrarily  close  to  1  for  sufficiently  small  5,  and 
because  of  (53),  for  any  positive  6 

ess  sup  \h{C)\  =  ess  sup \h{^)\ . 

|$|<5  ieR 

Assume  tp  fails  to  satisfy  (2).  To  see  that  this  impHes  that  A  also  fails  to 
satisfy  (2)  suppose,  on  the  contrary,  that  it  does  satisfy  (2).  Then 


(55) 


OO 

y:  |A(f-2rfr  =  i 

fe— —  OO 


and,  by  virtue  of  (54), 

(56) 

Thus,  we  may  write 

(57) 


ess  sup  \h{^)\  <  1 . 
ieR 

OO 

Y,  \H(  -  2Tk)\^  <  i{() 

k=—oo 
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where  $  is  the  trigonometric  polynomial 

OO 

k=—<x> 

Since  (p  fails  to  satisfy  (2)  the  polynomial  $(^)  is  arbitrarily  small  on  a  set 
of  positive  measure  which,  in  view  of  (57),  contradicts  (55). 

To  see  the  second  assertion  of  the  theorem  assume  that  A  satisfies  (2). 
Then  by  virtue  of  (55)  and  (56),  we  may  write 

(58)  K0i^(0r+  E  |^(^-27r/;)p>l. 

fce2\{o} 

Now,  since  $(^)  <  1  and  ^(0)  =  1,  we  may  conclude  that  for  any  e  satisfying 
0  <  e  <  1/2 

(59)  m)f  >1-6 

whenever  ^  is  in  a  sufficiently  small  neighborhood  of  0.  Since 

ess  inf  |h(£)l  =  ess  inf  |/i(£)| . 

1^|<«  ieR' 

for  any  positive  S,  (59)  allows  us  to  conclude  that 

ess  inf  |/i(OI  >  1. 

The  last  inequality  together  with  (56)  imply  that 

(60)  m)\  =  1 

almost  everywhere.  Conversely,  if  h  satisfies  (60)  and  $(^)  =  1  then  (55) 
follows  from  (52). 

2.7 

Proposition  3  is  an  easy  but  tedious  consequence  of  the  characterization  of 
zero  sets  of  those  polynomials  5'(^)  which  fail  to  generate  a  <p  which  satisfies 
(2).  These  sets  axe  characterized  in  [3,  5].  The  details  will  appear  elsewhere. 


3  Remarks  and  acknowledgements 

Section  1.1  was  motivated  by  Lawton’s  result  [6]  which  asserts  that  if  {sfc} 
is  a  finite  sequence  satisfying  (3)  and  (4)  and  if  y>  is  the  function  defined  by 
(5)  then  the  collection  of  functions  {2*'^^V’(2*®  —  j)}k,j&z  i  where 

(61)  V’(®)  =  -  k) 
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is  always  a  tight  frame.  See  also  [4,  page  178].  The  material  ia  Section  1.1 
impHes  that  there  are  also  orthonormal  wavelet  bases  associated  with  {s^} 
which  can  be  derived  from  the  scaling  function  (fa  in  the  usual  way. 

The  questions  addressed  in  this  paper  arose  in  lectures  WRM  presented 
during  the  1993  spring  semester  at  the  University  of  Connecticut.  Apparently 
the  question  raised  in  Subsection  1.4  was  first  publicly  raised  by  D.  PoUen, 
see  [3]. 

We  wish  to  thank  I.  Daubechies,  who  kindly  suggested  that  pathological 
scaling  sequences  which  axe  not  dilates  might  exist,  and  K.  Grochenig,  who 
brought  [3]  and  [5]  to  our  attention. 

Both  authors  were  partially  supported  by  DAPRA  Grant  AFOSR-90- 
0311. 
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